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Abstract 

We consider the problem of maintaining the Euclidean Delaunay triangulation DT of a set P of n 
moving points in the plane, along algebraic tranjectories of constant description complexity. Since the 
best known upper bound on the number of topological changes in the full Delaunay triangulation is only 
nearly cubic, we seek to maintain a suitable portion of the diagram that is less volatile yet retains many 
useful properties of the full triangulation. We introduce the notion of a stable Delaunay graph, which is 
a dynamic subgraph of the Delaunay triangulation. The stable Delaunay graph (a) is easy to define, (b) 
experiences only a nearly quadratic number of discrete changes, (c) is robust under small changes of the 
norm, and (d) possesses certain useful properties for further applications. 

The stable Delaunay graph (SDG in short) is defined in terms of a parameter a > 0, and consists of 
Delaunay edges pq for which the (equal) angles at which p and q see the corresponding Voronoi edge e pq 
are at least a. We show that (i) SDG always contains at least roughly one third of the Delaunay edges at 
any fixed time; (ii) it contains the /3-skeleton of P, for /3 = 1 + il(a 2 ); (iii) it is stable, in the sense that 
its edges survive for long periods of time, as long as the orientations of the segments connecting (nearby) 
points of P do not change by much; and (iv) stable Delaunay edges remain stable (with an appropriate 
redefinition of stability) if we replace the Euclidean norm by any sufficiently close norm. 

In particular, if we approximate the Euclidean norm by a polygonal norm (with a regular fc-gon as 
its unit ball, with k = 0(l/a)), we can define and keep track of a Euclidean SDG by maintaining the 
full Delaunay triangulation of P under the polygonal norm (which is trivial to do, and which is known 
to involve only a nearly quadratic number of discrete changes). 

We describe two kinetic data structures for maintaining SDG when the points of P move along 
pseudo-algebraic trajectories of constant description complexity. The first uses the polygonal norm ap- 
proximation noted above, and the second is slightly more involved, but significantly reduces the depen- 
dence of its performance on a. Both structures use O* (n) storage and process O* (n 2 ) events during the 
motion, each in 0*(1) time. (Here the 0*(-) notation hides multiplicative factors which are polynomial 
in 1/a and polylogarithmic in n.) 



*A preliminary version of this paper appeared in Proc. 26th Annual Symposium on Computational Geometry, 2010, pp. 127- 
136. 
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1 Introduction 



Delaunay triangulations and Voronoi diagrams. Let P be a (finite) set of points in M?. Let VD(P) and 
DT(P) denote the Voronoi diagram and Delaunay triangulation of P, respectively. For a point p G P, let 
Vor(p) denote the Voronoi cell of p. The Delaunay triangulation DT = DT(P) consists of all triangles 
whose circumcircles do not contain points of P in their interior. Its edges form the Delaunay graph, which 
is the straight-edge dual graph of the Voronoi diagram of P. That is, pq is an edge of the Delaunay graph if 
and only if Vor(p) and Vor (q) share an edge, which we denote by e pq . This is equivalent to the existence of 
a circle passing through p and q that does not contain any point of P in its interior — any circle centered at 
a point on e pq and passing through p and q is such a circle. Delaunay triangulations and Voronoi diagrams 
are fundamental to much of computational geometry and its applications. See (5j [TQ for a survey and a 
textbook on these structures. 

In many applications of Delaunay/Voronoi methods (e.g., mesh generation and kinetic collision de- 
tection) the points are moving continuously, so these diagrams need to be efficiently updated as motion 
occurs. Even though the motion of the nodes is continuous, the combinatorial and topological structure 
of the Voronoi and Delaunay diagrams change only at discrete times when certain critical events occur. 
Their evolution under motion can be studied within the framework of kinetic data structures (KDS in short) 
of Basch et al. [El El, a general methodology for designing efficient algorithms for maintaining such 
combinatorial attributes of mobile data. 

For the purpose of kinetic maintenance, Delaunay triangulations are nice structures, because, as men- 
tioned above, they admit local certifications associated with individual triangles. This makes it simple to 
maintain DT under point motion: an update is necessary only when one of these empty circumcircle con- 
ditions fails — this corresponds to cocircularities of certain subsets of four pointsQ Whenever such an event 
happens, a single edge flip easily restores Delaunayhood. Estimating the number of such events, however, 
has been elusive — the problem of bounding the number of combinatorial changes in DT for points moving 
along semi-algebraic trajectories of constant description complexity has been in the computational geometiy 
lore for many years; see ifTOl . 

Let n be the number of moving points in P. We assume that each point moves along an algebraic 
trajectory of fixed degree or, more generally, along pseudo- algebraic trajectory of constant description com- 
plexity (see Section|2]for a more formal definition). Guibas et al. |[T4l showed a roughly cubic upper bound 
of 0(n 2 X s (n)) on the number of discrete (also known as topological) changes in DT, where A s (n) is the 
maximum length of an (n, s)-Davenport-Schinzel sequence ll22l . and s is a constant depending on the mo- 
tions of the points. A substantial gap exists between this upper bound and the best known quadratic lower 
bound (22l. 

It is thus desirable to find approaches for maintaining a substantial portion of DT that provably expe- 
riences only a nearly quadratic number of discrete changes, that is reasonably easy to define and maintain, 
and that retains useful properties for further applications. 

Polygonal distance functions. If the "unit ball" of our underlying norm is polygonal then things improve 
considerably. In more detail, let Q be a convex polygon with a constant number, k, of edges. It induces a 
convex distance function 

dq(x, y) = min{A | y € x + XQ}; 

dQ is a metric if Q is centrally symmetric with respect to the origin. 

1 We assume that the motion of the points is sufficiently generic, so that no more than four points can become cocircular at any 
given time. 
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We can define the Q-Voronoi diagram of a set P of points in the plane in the usual way, as the par- 
titioning of the plane into Voronoi cells, so that the cell Vor°(p) of a point p is {x € M 2 \ cIq(x,p) = 
minp/gp cIq(x,p')}. Assuming that the points of P are in general position with respect to Q, these cells are 
nonempty, have pairwise disjoint interiors, and cover the plane. 

As in the Euclidean case, the Q-Voronoi diagram of P has its dual representation, which we refer to 
as the Q-Delaunay triangulation DT°(P) = DT°. A triple of points in P define a triangle in DT° if and 
only if they lie on the boundary of some homothetic copy of Q that does not contain any point of P in 
its interior. Assuming that P is in general position, these Q-Delaunay triangles form a triangulation of a 
certain simply-connected polygonal region that is contained in the convex hull of P. Unlike the Euclidean 
case, it does not always coincide with the convex hull (see Figures [5] and [23] for examples). See Chew and 
Drysdale (H and Leven and Shark |[T9l for analysis of Voronoi and Delaunay diagrams of this kind. 

For kinetic maintenance, polygonal Delaunay triangulations are "better" than Euclidean Delaunay tri- 
angulations because, as shown by Chew Q, when the points of P move (in the algebraic sense assumed 
above), the number of topological changes in the Q-Delaunay triangulation is only nearly quadratic in n. 
One of the major observations in this paper is that the stable portions of the Euclidean Delaunay triangula- 
tion and the Q-Delaunay triangulation are closely related. 



Stable Delaunay edges. We introduce the notion of a-stable Delaunay edges, for a fixed parameter a > 0, 
defined as follows. Let pq be a Delaunay edge under the Euclidean norm, and let Apqr + and Apqr~ be 
the two Delaunay triangles incident to pq. Then pq is called a-stable if its opposite angles in these triangles 
satisfy Zpr + q + Zpr~q < tt — a. (The case where pq lies on the convex hull of P is treated as if one of 
r + , r~ lies at infinity, so that the corresponding angle Zpr + q or Zpr~q is equal to 0.) An equivalent and 
more useful definition, in terms of the dual Voronoi diagram, is that pq is a-stable if the equal angles at 
which p and q see their common Voronoi edge e pq are at least a. See Figured] 




Figure 1: The points p and q see their common Voronoi edge ab at (equal) angles (3. This is equivalent to 
the angle condition x + y = n - f) for the two adjacent Delaunay triangles. 



A justification for calling such edges stable lies in the following observation: If a Delaunay edge pq 
is a-stable then it remains in DT during any continuous motion of P for which every angle Z.prq, for 
rGP\{p, q}, changes by at most a/2. This is clear because at the time pq is a-stable we have Zpr + q + 
/Lpr~q < tt — a for any pair of points r + , r~ lying on opposite sides of the line £ supporting pq, so, if 
each of these angles change by at most a/2 we still have Zpr + q + Zpr~q < tt, which is easily seen to 
imply that pq remains an edge of DT. (This argument also covers the cases when a point r crosses t from 
side to side: Since each point, on either side of I, sees pq at an angle of < tt — a, it follows that no point 
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can cross pq itself - the angle has to increase from tt — a to it. Any other crossing of £ by a point r causes 
Zprq to decrease to 0, and even if it increases to a/2 on the other side of I, pq is still an edge of DT, as is 
easily checked.) Hence, as long as the "small angle change" condition holds, stable Delaunay edges remain 
a "long time" in the triangulation. 

Informally speaking, the non-stable edges pq of DT are those for p and q are almost cocircular with 
their two common Delaunay neighbors r + , r~, and hence are more likely to get flipped "soon". 



Overview of our results. Let a > be a fixed parameter. In this paper we show how to maintain a 
subgraph of the full Delaunay triangulation DT, which we call a (ca, a)-stable Delaunay graph (SDG in 
short), so that (i) every edge of SDG is a-stable, and (ii) every ca-stable edge of DT belongs to SDG, where 
c > 1 is some (small) absolute constant. Note that SDG is not uniquely defined, even when c is fixed. 

In Section [2 we introduce several useful definitions and show that the number of discrete changes in 
the SDGs that we consider is nearly quadratic. What this analysis also implies is that if the true bound for 
kinetic changes in a Delaunay triangulation is really close to cubic, then the overhelming majority of these 
changes involve edges which never become stable and just flicker in and out of the diagram by cocircularity 
with their two Delaunay neighbors. 

In Sections [3] and 0] we show that SDG can be maintained by a kinetic data structure that uses only 
near-linear storage (in the terminology of ED, it is compact), encounters only a nearly quadratic number of 
critical events (it is efficient), and processes each event in polylogarithmic time (it is responsive). For the 
second data structure, described in Section @] can be slightly modified to ensure that each point appears at 
any time in only polylogarithmically many places in the structure (it then becomes local). 

The scheme described in Section|3]is based on a useful and interesting "equivalence" connection between 
the (Euclidean) SDG and a suitably defined "stable" version of the Delaunay triangulation of P under the 
"polygonal" norm whose unit ball Q is a regular fe-gon, for A: = 0(l/a). As noted above, Voronoi and 
Delaunay structures under polygonal norms are particularly favorable for kinetic maintenance because of 
Chew's result @, showing that the number of topological changes in DT°(P) is (9*(n 2 /c 4 ); here the 0*(-) 
notation hides a factor that depends sub-polynomially on both n and k. In other words, the scheme simply 
maintains the "polygonal" diagram DT°(P) in its entirety, and selects from it those edges that are also stable 
edges of the Euclidean diagram DT. 

The major disadvantage of the solution in Section |3]is the rather high (proportional to 6(l/a 4 )) depen- 
dence on l/a(~ k) of the bound on the number of topological changes. We do not know whether the upper 
bound 0*(n 2 k ) on the number of topological changes in DT°(P) is nearly tight (in its dependence on k). 
To remedy this, we present in Section |4]an alternative scheme for maintaining stable (Euclidean) Delaunay 
edges. The scheme is reminiscent of the kinetic schemes used in [[Q for maintaining closest pairs and nearest 
neighbors. It extracts 0*{n) pairs of points of P that are candidates for being stable Delaunay edges. Each 
point p € P then runs 0(1/ a) kinetic and dynamic tournaments involving the other points in its candidate 
pairs. Roughly, these tournaments correspond to shooting 0(1/ a) rays from P in fixed directions and find- 
ing along each ray the nearest point equally distant from p and from some other candidate point q. We show 
that pq is a stable Delaunay edge if and only if q wins many (at least some constant number of) consecutive 
tournaments of p (or p wins many consecutive tournaments of q). A careful analysis shows that the number 
of events that this scheme processes (and the overall processing time) is only 0*(n 2 /a 2 ). 

Section [5] establishes several useful properties of stable Delaunay graphs. In particular, we show that at 



any given time the stable subgraph contains at least 
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n Delaunay edges, i.e., at least about 



2(tt/q-2) 

one third of the maximum possible number of edges. In addition, we show that at any moment the SDG 
contains the closest pair, the so-called [3-skeleton of P, for [3 = 1 + £l(a 2 ) (see BETH), an d the crust of a 
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sufficiently densely sampled point set along a smooth curve (see GHU). We also extend the connection in 
Section [3] to arbitrary distance functions dq whose unit ball Q is sufficiently close (in the Hausdorff sense) 
to the Euclidean one (i.e., the unit disk). 



2 Preliminaries 

Stable edges in Voronoi diagrams. Let {no, • • • ,Uk-i} C S 1 be a set of k = 0(1/ a) equally spaced 
directions in M 2 . For concreteness take u- L = (cos(2m/k), — sm(2m/k)), < i < k (so our directions Uj 
go clockwise as i increases) H For a point p € P and a unit vector u let u[p] denote the ray {p + \u | A > 0} 
that emanates from p in direction u. For a pair of points p,q £ P let b pq denote the perpendicular bisector 
of p and q. If b pq intersects Uj [p] , then the expression 

ii ii2 

is the distance between p and the intersection point of b pq with ui [p] . If b pq does not intersect Uj [p] we define 
PilPiOl = °°- The point q minimizes <fii\p,q'], among all points q' for which b pq > intersects Ui\p], if and 
only if the intersection between b pq and Ui[p] lies on the Voronoi edge e pq . We call q the neighbor of p in 
direction m, and denote it by Ni(p); see Figure |2] 

The (angular) extent of a Voronoi edge e pq of two points p, q G P is the angle at which it is seen from 
either p or q (these two angles are equal). For a given angle a < ir, e pq is called a-long (resp., a-short) if 
the extent of e pq is at least (resp., smaller than) a. We also say that pq E DT(P) is a-long (resp., a-short) if 
e pq is a-long (resp., a-short). As noted in the introduction, these notions can also be denned (equivalently) 
in terms of the angles in the Delaunay triangulation: A Delaunay edge pq, which is not a hull edge, is a-long 
if and only if Zpr + q + Zpr~q < it — a, where Apr + q and Apr^q are the two Delaunay triangles incident 
to pq. See Figure [TJ hull edges are handled similarly, as discussed in the introduction. 

Given parameters a' > a > 0, we seek to construct (and maintain under motion) an (a' , a)-stable 
Delaunay graph (or stable Delaunay graph, for brevity, which we further abbreviate as SDG) of P, which 
is any subgraph G of DT(P) with the following properties: 

(51) Every a'-long edge of DT(P) is an edge of G. 

(52) Every edge of G is an a-long edge of DT(P). 

An (a', a)-stable Delaunay graph need not be unique. In what follows, a 1 will always be some fixed (and 
reasonably small) multiple of a. 



Kinetic tournaments. Kinetic tournaments were first studied by Basch et al. O, for kinetically main- 
taining the lowest point in a set P of n points moving on some vertical line, say the y-axis, so that their 
trajectories are algebraic of bounded degree, as above. These tournaments are a key ingredient in the data 
structures that we will develop for maintaining stable Delaunay graphs. Such a tournament is represented 
and maintained using the following variant of a heap. Let T be a minimum-height balanced binary tree, 
with the points stored at its leaves (in an arbitrary order). For an internal node v G T, let P v denote the 
set of points stored in the subtree rooted at v. At any specific time t, each internal node v stores the lowest 

2 The index arithmetic is modulo k, i.e., m = Ui+k- 
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point among the points in P v at time t, which is called the winner at v. The winner at the root is the desired 
overall lowest point of P. 

To maintain T we associate a certificate with each internal node v, which asserts which of the two 
winners, at the left child and at the right child of v, is the winner at v. This certificate remains valid as long 
as (i) the winners at the children of v do not change, and (ii) the order along the y-axis between these two 
"sub-winners" does not change. The actual certificate caters only to the second condition; the first will be 
taken care of recursively. Each certificate has an associated failure time, which is the next time when these 
two winners switch their order along the y-axis. We store all certificates in another heap, using the failure 
times as keysJl This heap of certificates is called the event queue. 

Processing an event is simple. When the two sub-winners p, q at some node v change their order, we 
compute the new failure time of the certificate at v (the first future time when p and q meet again), update 
the event queue accordingly, and propagate the new winner, say p, up the tree, revising the certificates at the 
ancestors of v , if needed. 

If we assume that the trajectories of each pair of points intersect at most r times then the overall number 
of changes of winners, and therefore also the overall number of events, is at most ^ v \P(v)\/3 r (\P(v)\) = 
0(nf3 r (n) log n). Here /3 r (n) = X r (n)/n, and A r (n) is the maximum length of a Davenport-Schinzel 
sequence of order r on n symbols; see ll22l . 

This is larger by a logarithmic factor than the maximum possible number of times the lowest point along 
the y-axis can indeed change, since this latter number is bounded by the complexity of the lower envelope 
of the trajectories of the points in P (which, as noted above, records the changes in the winner at the root of 
T). 

Agarwal et al. HI show how to make such a tournament also dynamic, supporting insertions and dele- 
tions of points. They replace the balanced binary tree T by a weight-balanced (BB(a)) tree ||2T1 (and see 
also [20]). This allows us to insert a new point anywhere we wish in T, and to delete any point from T, 
in O(logn) time. Each such insertion or deletion may change O(logn) certificates, along the correspond- 
ing search path, and therefore updating the event queue takes 0(log 2 n) time, including the time for the 
structural updates of (rotations in) T; here n denotes the actual number of points in T, at the step where we 
perform the insertion or deletion. The analysis of [JT] is summarized in Theorem 12.11 

Theorem 2.1 (Agarwal et al. HI). A sequence of m insertions and deletions into a kinetic tournament, 
whose maximum size at any time is n (assuming m > n), when implemented as a weight-balanced tree in 
the manner described above, generates at most 0{mf3 r+ 2{n) logn) events, with a total processing cost of 
0(m/3 r _)_2( n ) log 2 n). Here r is the maximum number of times a pair of points intersect, and P r +2( n ) = 
\ T+ 2{n)/n. Processing an update or a tournament event takes 0(log 2 n) worst-case time. A dynamic 
kinetic tournament on n elements can be constructed in 0(n) time. 

Remarks: (1) Theorem 12. 1 1 subsumes the static case too, by inserting all the elements "at the beginning of 
time", and then tracing the kinetic changes. 

(ii) Note that the amortized cost of an update or of processing a tournament event is only O(logn) (as 
opposed to the 0(log 2 n) worst-case cost). 

Maintenance of an SDG. LetP = {pi, . . . , p n } be a set of points moving in M. 2 . Letpi(t) = (xi(t) , yi(t)) 
denote the position of pi at time t. We call the motion of P algebraic if each Xi (t) , yi (i) is a polynomial 
function of t, and the degree of motion of P is the maximum degree of these polynomials. Throughout this 
paper we assume that the motion of P is algebraic and that its degree is bounded by a constant. In this 

3 Any "standard" heap that supports insert, delete, and deletemin in 0(log n) time is good for our purpose. 



5 



subsection we present a simple technique for maintaining a (2a, a)-stable Delaunay graph. Unfortunately 
this algorithm requires quadratic space. It is based on the following easy observation (see Figure 0), where 
k is an integer, and the unit vectors (directions) uo, . . . , Uk-\ are as defined earlier. 

Lemma 2.2. Let a = 2n/k. ( i) If the extent of e pq is larger than 2a then there are two consecutive directions 
Ui, Ui + \, such that q is the neighbor ofp in directions U{ and Uj+i. 

(ii) If there are two consecutive directions u.j, Uj+i, such that q is the neighbor ofp in both directions U{ and 
Ui+i, then the extent of e pq is at least a. 




Figure 2: q is the neighbor of p in the directions m and u i+1 , so the Voronoi edge e pq is a-long. 



The algorithm maintains Delaunay edges pq such that there are two consecutive directions Uj and Uj+i 
along which q is the neighbor of p. For each point p and direction m we get a set of at most n — 1 piecewise 
continuous functions of time, tp%\p, q], one for each point q ^ p, as defined in (Q]). (Recall that ifi[p, q] = oo 
when m [p] does not intersect b pq .) By assumption on the motion of P, for each p and q, the domain in which 
ifi [p, q] (t) is defined consists of a constant number of intervals. 

For each point p, and ray m [p], consider each function ifi [p, q] as the trajectory of a point moving along 
the ray and corresponding to q. The algorithm maintains these points in a dynamic and kinetic tournament 
Ki{p) (see Theorem 12. Il l that keeps track of the minimum of {<Pi\p,q)(t)} q7 L p over time. For each pair 
of points p and q such that q wins in two consecutive tournaments, Ki{p) and Ki + \(p), of p, it keeps the 
edge pq in the stable Delaunay graph. It is trivial to update this graph as a by-product of the updates of the 
various tournaments. The analysis of this data structure is straightforward using Theorem 12.11 and yields 
the following result. 

Theorem 2.3. Let P be a set of n moving points in M 2 under algebraic motion of bounded degree, let k 
be an integer, and let a = 2n/k. A (2a, a)-stable Delaunay graph of P can be maintained using 0(kn 2 ) 
storage and processing 0(kn 2 /3 r+ 2{n) logn) events, for a total cost of 0(kn 2 f3 r+ 2(n) log 2 n) time. The 
processing of each event takes 0(log 2 n) worst-case time. Here r is a constant that depends on the degree 
of motion of P. 

Later on, in Section HJ we will revise this approach and reduce the storage to nearly linear, by letting 
only a small number of points to participate in each tournament. The filtering procedure for the points makes 
the improved solution somewhat more involved. 
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3 An SDG Based on Polygonal Voronoi Diagrams 



Let Q = Qk be a regular fc-gon for some even k = 2s, circumscribed by the unit disk, and let a = 7r/s (this 
is the angle at which the center of Q sees an edge). Let VD°(P) and DT°(P) denote the Q-Voronoi diagram 
and the dual Q-Delaunay triangulation of P, respectively. In this section we show that the set of edges 
of VD°(P) with sufficiently many breakpoints (see below for details) form a /3') -stable (Euclidean) 
Delaunay graph for appropriate multiples j3,j3' of a. Thus, by kinetically maintaining VD°(P) (in its 
entirety), we shall get "for free" a KDS for keeping track of a stable portion of the Euclidean DT. 

3.1 Properties of VD°(P) 

We first review the properties of the (stationary) VD°(P) and DT°(P). Then we consider the kinetic version 
of these diagrams, as the points of P move, and review Chew's proof Q that the number of topological 
changes in these diagrams, over time, is only nearly quadratic in n. Finally, we present a straightforward 
kinetic data structure for maintaining DT°(P) under motion that uses linear storage, and that processes a 
nearly quadratic number of events, each in 0(log n) time. Although later on we will take Q to be a regular 
A;-gon, the analysis in this subsection is more general, and we only assume here that Q is an arbitrary convex 
fc-gon, lying in general position with respect to P. 

Stationary Q-diagrams. The bisector b^ q between two points p and q, with respect to <1q (■,•), is the locus 
of all placements of the center of any homothetic copy Q' of Q that touches p and q. Q' can be classified 
according to the pair of its edges, e\ and e2, that touch p and q, respectively. If we slide Q' so that its 
center moves along b° pq (and its size expands or shrinks to keep it touching p and q), and the contact edges, 
ei and e2, remain fixed, the center traces a straight segment. The bisector is a concatenation of 0(k) such 
segments. They meet at breakpoints, which are placements of the center of a copy Q' that touches p and q 
and one of the contact points is a vertex of Q; see Figure [3] We call such a placement a corner contact at 
the appropriate point. Note that a corner contact where some vertex w of (a copy Q' of) Q touches p has 
the property that the center of Q' lies on the fixed ray emanating from p and parallel to the directed segment 
from iu to the center of Q. 




Figure 3: Each breakpoint on b° pq corresponds to a corner contact of Q at one of the points p, q, so that dQ 
also touches the other point. 

A useful property of bisectors and Delaunay edges, in the special case where Q is a regular fc-gon, which 
will be used in the next subsection, is that the breakpoints along a bisector b^ q alternate between corner 
contacts at p and corner contacts at q. Indeed, assuming general position, each point w £ dQ determines 
a unique placement of Q where it touches p at w and also touches q, as is easily checked. A symmetric 
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property holds when we interchange p and q. Hence, as we slide the center of Q along the bisector b pq , 
the points of contact of dQ with p and q vary continuously and monotonically along dQ. Consider two 
consecutive corner contacts, Q' , Q" , of Q at p along tit „, and suppose to the contrary that the portion of b pq 
between them is a straight segment, meaning that, within this portion, dQ touches each of p, q at a fixed 
edge. Since the center of Q moves along the angle bisector of the lines supporting these edges (a property 
that is easily seen to hold for regular fc-gons), it is easy to see that the distance between the two contact 
points of p, at the beginning and the end of this sliding, and the distance between the two contact points of q 
(measured, say, on the boundary of the standard placement of Q) are equal. However, this distance for p is 
the length of a full edge of dQ, because the motion starts and ends with p touching a vertex, and therefore 
the same holds for q, which is impossible (unless q also starts and ends at a vertex, which contradicts our 
general position assumption). 

Another well known property of Q-bisectors and Voronoi edges, for arbitrary convex polygons in general 
position with respect to P, is that two bisectors b pqi , b pq2 , can intersect at most once (again, assuming 
general position), so every Q-Voronoi edge e pq is connected. Equivalently, this asserts that there exists at 
most one homothetic placement of Q at which it touches p, q\, and q^. Indeed, since homothetic placements 
of Q behave like pseudo-disks (see, e.g., ifTTI ). the boundaries of two distinct homothetic placements of Q 
intersect in at most two points, or, in degenerate position, in at most two connected segments. Clearly, in the 
former case the boundaries cannot both contain p, q\, and q2 , and this also holds in the latter case because 
of our general position assumption. 

Consider next an edge pq of DT°(P). Its dual Voronoi edge e pq is a portion of the bisector b pq , and 
consists of those center placements along b pq for which the corresponding copy Q' has an empty interior 
(i.e., its interior is disjoint from P). Following the notation of Chew Q, we call pq a corner edge if e p 
contains a breakpoint (i.e., a placement with a corner contact); otherwise it is a non-corner edge, and is 
therefore a straight segment. 

Kinetic Q _ diagrams. Consider next what happens to VD°(P) and DT°(P) as the points of P move 
continuously with time. In this case VD°(P) changes continuously, but undergoes topological changes at 
certain critical times, called events. There are two kinds of events: 

(i) Flip Event. A Voronoi edge e pq shrinks to a point, disappears, and is "flipped" into a newly emerging 
Voronoi edge e plql . 

(ii) Corner Event. An endpoint of some Voronoi edge e pq becomes a breakpoint (a comer placement). 
Immediately after this time e pq either gains a new straight segment, or loses such a segment, that it had 
before the event. 

Some comments are in order: 

(a) A flip event occurs when the four points p, q,p' , q' become "cocircular": there is an empty homothetic 
copy Q' of Q that touches all four points. 

(b) Only non-corner edges can participate in a flip event, as both the vanishing edge e pq and the newly 
emerging edge e p , , do not have breakpoints near the event. 

(c) A flip event entails a discrete change in the Delaunay triangulation, whereas a corner event does not. 
Still, for algorithmic puiposes, we will keep track of both kinds of events. 

We first bound the number of corner events. 

Lemma 3.1. Let P be a set of n points in M 2 under algebraic motion of bounded degree, and let Q be a 
convex k-gon. The number of corner events in DT°(P) is 0(k 2 nX r (n)), where r is a constant that depends 
on the degree of motion of P. 
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Proof. Fix a point p and a vertex w of Q, and consider all the corner events in which w touches p. As noted 
above, at any such event the center c of Q lies on a ray 7 emanating from p at a fixed direction. (Since 
p is moving, 7 is a moving ray, but its orientation remains fixed.) For each other point q € P \ {p}, let 
cpZ [p, q] denote the distance, at time t, from p along 7 to the center of a copy of Q that touches p (at w) and 
q. The value min g <£>*[p, (?](t) represents the intersection of d Vor°(p) with 7 at time where Vor°(p) is the 
Voronoi cell of p in VD°(P). The point q that attains the minimum defines the Voronoi edge ei (or vertex 
if the minimum is attained by more than one point q) of Vor (p) that 7 intersects. 

In other words, we have a collection of re — 1 partially defined functions ip°\p, q], and the breakpoints 
of their lower envelope represent the comer events that involve the contact of w with p. By our assumption 
on the motion of P, each function ip°\p, q] is piecewise algebraic, with 0(k) pieces. Each piece encodes 
a continuous contact of q with a specific edge of Q', and has constant description complexity. Hence (see, 
e.g., E2l Corollary 1.6]) the complexity of the envelope is at most 0(kX r (n)), for an appropriate constant 
r. Repeating the analysis for each point p and each vertex w of Q, the lemma follows. □ 

Consider next flip events. As noted, each flip event involves a placement of an empty homothetic copy 
Q' of Q that touches simultaneously four points pi,P2,P3,P4 of P, in this counterclockwise order along 
dQ', so that the Voronoi edge e pip3 , which is a non-corner edge before the event, shrinks to a point and is 
replaced by the non-corner edge e p2P4 - Let denote the edge of Q' that touches pi, for i = 1, 2, 3, 4. 

We fix the quadruple of edges ei,e2,e3,e4, bound the number of flip events involving a quadruple 
contact with these edges, and sum the bound over all 0(k 4 ) choices of four edges of Q. For a fixed quadruple 
of edges ei,e2,e3,e4, we replace Q by the convex hull Qo of these edges, and note that any flip event 
involving these edges is also a flip event for Qo- We therefore restrict our attention to Qq, which is a convex 
fco-gon, for some ko < 8. 

We note that if (p, q) is a Delaunay edge representing a contact of some homothetic copy Q' of Qo 
where p and q touch two adjacent edges of Q' , then (p, q) must be a corner edge — shrinking Q' Q towards 
the vertex common to the two edges, while it continues to touch p and q, will keep it empty, and eventually 
reach a placement where either p or q touches a corner of Q' Q . The same (and actually simpler) argument 
applies to the case when p and q touch the same edg^l of Qo- 




Figure 4: Left: The edge ef 3 in the diagram VD°(P) before disappearing. The endpoint c1 23 (resp., c* 43 ) 
of e* 3 corresponds to the homothetic copy of Qo whose edges e%, e 2 , e 3 (resp., ei, e 4 , e 3 ) are incident to the 
respective vertices pi,p2,P3 (resp., pi,pa,P3). Right: The tree of non-corner edges. 

Consider the situation just before the critical event takes place, as depicted in Figure [4] (left). The 
Voronoi edge ep ip3 (to simplify the notation, we write this edge as ef 3 , and similarly for the other edges and 

4 In general position this does not occur, but it can happen at discrete time instances during the motion, 
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vertices in this analysis) is delimited by two Voronoi vertices, one, c\ 2 ^, being the center of a copy of Qq 
which touches pi,P2,P3 at the respective edges e-y, e^, and the other, cf 43 , being the center of a copy of 
Qo which touches pi,P4,,P3 at the respective edges &\^e^&z. Consider the two other Voronoi edges e\ 2 and 
e 23 adjacent to c 423 , and the two Voronoi edges ef 4 and e 43 adjacent to cf 43 . Among them, consider only 
those which are non-corner edges; assume for simplicity that they all are. For specificity, consider the edge 
e\ 2 . As we move the center of Qq along that edge away from c^ 23 , Qo loses the contact with p 3 ; it shrinks 
on the side of P1P2 which contains p 3 (and p 4 , already away from Qq), and expands on the other side. Since 
this is a non-comer edge, its other endpoint is a placement where the (artificial) edge e\2 of Qo between e\ 
and e 2 touches another point p§. Now, however, since e±2 is adjacent to both edges e±, 62, the new Voronoi 
edges ef 5 and e 25 are both corner edges. 

Repeating this analysis to each of the other three Voronoi edges adjacent to e* 3 , we get a tree of non- 
corner Voronoi edges, consisting of at most five edges, so that all the other Voronoi edges adjacent to its 
edges are corner edges. As long as no discrete change occurs at any of the surrounding corner edges, the 
tree can undergo only 0(1) discrete changes, because all its edges are defined by a total of O(l) points of 
P. When a corner edge undergoes a discrete change, this can affect only O(l) adjacent non-corner trees of 
the above kind. Hence, the number of changes in non-corner edges is proportional to the number of changes 
in corner edges, which, by Lemma 137X1 (applied to Qo) is 0(nX r (n)). Multiplying by the 0{k i ) choices of 
quadruples of edges of Q, we thus obtain: 

Theorem 3.2. Let P be a set of n moving points in M 2 under algebraic motion of bounded degree, and let 
Q be a convex k-gon. The number of topological changes in VD°(P) with respect to Q is 0(k 4 n\ r (n)), 
where r is a constant that depends on the degree of motion of P. 

Kinetic maintenance of VD°(P) and DT°(P). As already mentioned, it is a fairly trivial task to main- 
tain DT°(P) and VD°(P) kinetically, as the points of P move. All we need to do is to assert the correctness 
of the present triangulation by a collection of local certificates, one for each edge of the diagram, where the 
certificate of an edge asserts that the two homothetic placements Q~ , Q + of Q that circumscribe the two 
respective adjacent Q-Delaunay triangles Apqr~ , Apqr + , are such that Q~ does not contain r + and Q + 
does not contain r~. The failure time of this certificate is the first time (if one exists) at which p, q, r~, and 
r + become O-cocircular — they all lie on the boundary of a common homothetic copy of Q. Such an event 
corresponds to a flip event in DT°(P). If pq is an edge of the periphery of DT°(P), so that Apqr + exists 
but Apqr~ does not, then Q~ is a limiting wedge bounded by rays supporting two consecutive edges of (a 
copy of) Q, one passing through p and one through q (see Figure [5]). The failure time of the corresponding 
certificate is the first time (if any) at which r + also lies on the boundary of that wedge. 

We maintain the breakpoints using "sub-certificates", each of which asserts that say, touches each 
of p, q, r~ at respective specific edges (and similarly for Q + ). The failure time of this sub-certificate is the 
first failure time when one of p, q or r~ touches Q~ at a vertex. In this case we have a corner event — two of 
the adjacent Voronoi edges terminate at a corner placement. Note that the failure time of each sub-certificate 
can be computed in 0(1) time. Moreover, for a fixed collection of valid sub-certificates, the failure time of 
an initial certificate (asserting non-cocircularity) can also be computed in O(l) time (provided that it fails 
before the failures of the corresponding sub-certificates), because we know the four edges of Q~ involved 
in the contacts. 

We therefore maintain an event queue that stores and updates all the active failure times (there are only 
O(n) of them at any given time — the bound is independent of k, because they correspond to actual DT 
edges. When a sub-certificate fails we do not change DT°(i- > ), but only update the corresponding Voronoi 
edge, by adding or removing a segment and a breakpoint, and by replacing the sub-certificate by a new one; 
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Figure 5: If r does not exist then Q is a limiting wedge bounded by rays supporting two consecutive 
edges of (a copy of) Q. 



we also update the cocircularity certificate associated with the edge, because one of the contact edges has 
changed. When a cocircularity certificate fails we update DT°(P) and construct 0(1) new sub-certificates 
and certificates. Altogether, each update of the diagram takes 0(log n) time. We thus have 

Theorem 3.3. Let P be a set of n moving points in M? under algebraic motion of bounded degree, and let 
Q be a convex k-gon. DT°(P) and VD°(P) can be maintained using 0(n) storage and O(logn) update 
time, so that 0(k 4 nX r (n)) events are processed, where r is a constant that depends on the degree of motion 
of P. 

3.2 Stable Delaunay edges in DT°(P) 

We now restrict Q to be a regular k-gon. Let vq, . . . , Vk-i be the vertices of Q arranged in a clockwise 
direction, with vq the leftmost. We call a homothetic copy of Q whose vertex vj touches a point p, a Vj- 
placement of Q at p. Let Uj be the direction of the vector that connects Vj with the center of Q, for each 
< j < k (as in Section|2]>. See Figure |6](left). 

We follow the machinery in the proof of Lemma I3TI That is, for any pair p,q € P let </?|fai q] denote the 
distance from p to the point Uj [p] n bp q ; we put [p, q] = oo if Uj [p] does not intersect bp q . If tp? [p, q] < oo 
then the point b^ q n Uj[p] is the center of the Vj -placement Q' of Q at p that also touches q, and it is easy 
to see that there is a unique such point. The value ipj[q,p] is equal to the circumradius of Q' . See Figure [6] 
(middle). 

The neighbor Nj\p] of p in direction Uj is defined to be the point q € P\ {p} that minimizes tp? [p, q]. 
Clearly, for any p,q € P we have iVJ[p] = q if and only if there is an empty Vj -placement Q' of Q at p so 
that q touches one of its edges. 

Remark: Note that, in the Euclidean case, we have ipj [p, q] < oo if and only if the angle between pq and 
Uj\p] is at most tt/2. In contrast, fj\p, q] < oo if and only if the angle between pq and Uj[p] is at most 
7r/2 — it /k = it /2 — a/2. Moreover, we have <p>j\p, q] < Pj[p, q]. Therefore, </?*[p, q] < oo always implies 
ipj [p, q] < oo, but not vice versa; see Figure |6](right). Note also that in both the Euclidean and the polygonal 
cases, the respective quantities Nj [p] and [p] may be undefined. 

Lemma 3.4. Let p,q € P be a pair of points such that Nj(p) = qfor h > 3 consecutive indices, say 
< j < h — 1. Then for each of these indices, except possibly for the first and the last one, we also have 
Nf]p] = q. 

Proof. Let w\ (resp., W2) be the point at which the ray Uo[p] (resp., %_i[p]) hits the edge e pq in VD(P). 
(By assumption, both points exist.) Let D\ and D2 be the disks centered at w\ and W2, respectively, and 
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Figure 6: Left: Uj is the direction of the vector connecting vertex Vj to the center of Q. Middle: The function 
ipjlp, q] is equal to the radius of the circle that circumscribes the Uj-placement of Q at p that also touches q. 
Right: The case when <p$\p, q] = oo while <fj\p,q] < oo. In this case q must lie in one of the shaded wedges. 



touching p and q. By definition, neither of these disks contains a point of P in its interior. The angle between 
the tangents to D\ and D2 at p or at q (these angles are equal) is j3 = (h — l)a; see Figure [7J (left). 




Figure 7: Left: The angle between the tangents to D x and D 2 atp(or at q) is equal to Zw!pw 2 = /? = (h-l)a. 
Right: The line I' crosses D in a chord qq' which is fully contained in e'. 

Fix an arbitrary index 1 < j < h — 2, so nj [p] intersects e pg and forms an angle of at least a with 
each of pwi,pw 2 - Let Q' be the Uj-placement of Q at p that touches q. To see that such a placement exists, 
we note that, by the preceding remark, it suffices to show that the angle between pq and Uj [p] is at most 
7r/2 — a/2; that is, to rule out the case where q lies in one of the shaded wedges in Figure |6] (right). This 
case is indeed impossible, because then one of Uj-\[p},Uj + i[p\ would form an angle greater than tt/2 with 
pq, contradicting the assumption that both of these rays intersect the (Euclidean) b pq . 

We claim that Q' C D\ U D2. Establishing this property for every 1 < j < h— 2 will complete the proof 
of the lemma. Let e' be the edge of Q' passing through q. See Figure [7] (right). Let D be the disk whose 
center lies on Uj\p] and which passes through p and q, and let D + be the circumscribing disk of Q'. Since 
q € dD and is interior to D + , and since D and D + are centered on the same ray Uj[q] and pass through 
p, it follows that D C D + . The line £' containing e' crosses D in a chord qq' that is fully contained in e'. 
The angle between the tangent to D at q and the chord qq' is equal to the angle at which p sees qq'. This is 
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smaller than the angle at which p sees e', which in turn is equal to a/2. 

Arguing as in the analysis of D\ and D2, the tangent to D at q forms an angle of at least a with each 
of the tangents to Di, D2 at q, and hence e' forms an angle of at least a/2 with each of these tangents; see 
Figure [8](left). The line £' marks two chords gig, gg2 within the respective disks D±, D2. We claim that e' is 
fully contained in their union gig2- Indeed, the angle qipq is equal to the angle between I' and the tangent 
to D\ at q, so Zgipg > a/2. On the other hand, the angle at which p sees e' is a/2, which is smaller. This, 
and the symmetic argument involving D2, are easily seen to imply the claim. 




Figure 8: Left: The line £' forms an angle of at least a/2 with each of the tangents to D 1 ,D 2 at q. Right: 
The edge e' = aia 2 of Q' is fully contained in Di U D 2 . 



Now consider the circumscribing disk D + of Q'. Denote the endpoints of e' as a± and 02, where a± lies 
in gig and 02 lies in gg2. Since the ray pa 1 hits dD + before hitting D\, and the ray pq hits these circles in 
the reverse order, it follows that the second intersection of dD\ and dD + (other than p) must lie on a ray 
from p which lies between the rays pai,pq and thus crosses e'. See Figure [8] (right). Symmetrically, the 
second intersection point of dD2 and dD + also lies on a ray which crosses e' . 

It follows that the arc of dD + delimited by these intersections and containing p is fully contained in 
D\ U Z?2- Hence all the vertices of Q' (which lie on this arc) lie in D\ U D2. This, combined with the 
argument in the preceding paragraphs, is easily seen to imply that Q' C D\ U D2, so its interior does 
not contain points of P, which in turn implies that Nj\p] = q. As noted, this completes the proof of the 
lemma. □ 

Since Q-Voronoi edges are connected, Lemma l3T4l implies that is "long", in the sense that it contains 
at least h — 2 breakpoints that represent corner placements at p, interleaved (as promised in Section |3TTT > with 
at least h — 3 corner placements at q. This property is easily seen to hold also under the weaker assumptions 
that: (i) for the first and the last indices j = 0, h — 1, the point Nj[p] either is equal to q or is undefined, 
and (ii) for the rest of the indices j we have Nj\p] = q and ip? [p, q] < 00 (i.e., the vj -placement of Q at 
p that touches q exists). In this relaxed setting, it is now possible that any of the two points w\, W2 lies at 
infinity, in which case the corresponding disk D\ or D2 degenerates into a halfplane. This stronger version 
of Lemma l3.4l is used in the proof of the converse Lemma 1331 asserting that every edge e|„ in VD°(P) with 
sufficiently many breakpoints has a stable counterpart e pq in VD(P). 

Lemma 3.5. Let p,q € P be a pair of points such that Nj[p] = qfor at least three consecutive indices 
j € {0, . . . , k — 1}. Then for each of these indices, except possibly for the first and the last one, we have 
N j [p]=q. 

Proof. Again, assume with no loss of generality that Nj[p] = q for < j < h — 1, with h > 3. Suppose 
to the contrary that, for some 1 < j < h — 2, we have Nj [p] ^ q. Since [p] = q by assumption, we 
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have ipj \p, q] < </?|[p, q] < oo, so there exists r € P for which <£j [p, r] < ipj [p, q\. Assume with no loss of 
generality that r lies to the left of the line from p to q. In this case ipj_i\p, r] < ipj-\\p, q] < oo. Indeed, 
we have (i) Nj_i\p] = q by assumption, so <p*_i\p,q] < oo, and (ii) ipj_i\p,q] < (Pj_x\p,q\. Moreover, 
because r lies to the left of the line from p to q, the orientation of b pr lies counterclockwise to that of b pq , 
implying that ^_i[p, g] < oo. See Figure [9] Since Uj\p] hits b pr before hitting b pq , any ray emanating from 
p counterlockwise to uj[p] must do the same, so we have tpj_i\p, r] < ipj_i\p, q], as claimed. Similarly, 
we get that either <£>j_2[p, r] < tpj-2\p, q] < oo or ipj^ip, f] < ipj-2\p, q] = oo (where the latter can 
occur only for j = 1). Now applying (the extended version of) Lemma [3741 to the point set {p, q,r} and 
to the index set {j — 2,j — l,j}, we get that r] < (Pj_i\p, q\. But this contradicts the fact that 

Nf_ 1 \p] = q- □ 




Figure 9: Proof of Lemma l375l If Nj[p] ^ q because some r, lying to the left of the line from p to r, satisfies 

<Pj\p,r] < <pj\p,q]. Since ^_i[p,g] < ^|_i[p,g] < oo, we have ^_i[p,r] < ipj-^q}. 



Maintaining an SDG using VD°(P). Lemmas 13.41 and [3751 together imply that an SDG can be main- 
tained using the fairly straightforward kinetic algorithm for maintaining the whole VD°(P), provided by 
Theorem 13.31 We use VD°(P) to maintain the graph G on P, whose edges are all the pairs (p,q) € P x P 
such that p and q define an edge e pg in VD°(P) that contains at least seven breakpoints. As shown in 
Theorem 13.31 this can be done with 0(n) storage, O(logn) update time, and 0{k i n\ r {n)) updates (for an 
appropriate r). We claim that G is a (6a, a)-SDG in the Euclidean norm. 

Indeed, if two points p,q G P define a 6a-long edge e pq in VD(P) then this edge stabs at least six 
rays Uj\p] emanating from p, and at least six rays Uj[q] emanating from q. Thus, according to Lemma 13741 
VD°(P) contains the edge e pq with at least four breakpoints corresponding to corner placements of Q at 
p that touch q, and at least four breakpoints corresponding to comer placements of Q at q that touch p. 
Therefore, e pq contains at least 8 breakpoints, so (p, q) € G. 

For the second part, if p, q G P define an edge e pq in VD°(P) with at least 7 breakpoints then, by the 
interleaving property of breakpoints, we may assume, without loss of generality, that at least four of these 
breakpoints correspond to P-empty corner placements of Q at p that touch q. Thus, Lemma [3751 implies that 
VD(P) contains the edge e pq , and that this edge is hit by at least two consecutive rays Uj\p\. But then, as 
observed in Lemma 12721 the edge e pq is a-long in VD(P). We thus obtain the main result of this section. 

Theorem 3.6. Let P be a set of n moving points in M? under algebraic motion of bounded degree, and let 
a > be a parameter. A (6a, a)-stable Delaunay graph of P can be maintained by a KDS of linear size 
that processes 0(nX r (n)/a : ) events, where r is a constant that depends on the degree of motion of P, and 
that updates the SDG at each event in 0(log n) time. 
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4 An Improved Data Structure 



The data structure of Theorem 13 . 6 I requires 0(n) storage but the best bound we have on the number of events 
it may encounter is 0*(n 2 /a 4 ), which is much larger than the number of events encountered by the data 
structure of Theorem 12.3 1 (which, in terms of the dependence on a, is only 0*(n 2 /a)). In this section we 
present an alternative data structure that requires 0*(n/a 2 ) space and 0*(n 2 /a 2 ) overall processing time. 
The structure processes each event in 0*(l/a) time and is also local, in the sense that each point is stored 
at only 0*((l/a) 2 ) places in the structure. 

Notation. We use the directions u% and the associated quantities JVj [p] and tpi [p, q] defined in Section [2] 
We assume that k, the number of canonical directions, is even, and write, as in Section |2j k = 2s. We 
denote by Cj the cone (or wedge) with apex at the origin that is bounded by u, and Ui + \. Note that C, and 
Ci± s are antipodal. As before, for a vector it, we denote by u[x] the ray emanating from x in direction u. 
Similarly, for a cone C we denote by C[x] the translation of C that places its apex at x. Let 0</3<7r/2be 
an angle. For a direction u G S 1 and for two points p,q G P, we say that the edge e pq G VD(P) is (3-long 
around the ray u[q] if p is the Voronoi neighbor of q in all directions in the range [u — f3,u + 0\, i.e., for all 
v G [u — j3, u + 0\, the ray v[q] intersects e pq . The (5 -cone around u[q] is the cone whose apex is q and each 
of its bounding rays makes an angle of (3 with u[q\. 




Figure 10: Left: q is j-extremal for p. Right: q is strongly j-extremal for p. 



Definition (j'-extremal points), (i) Let p,q G P, let i be the index such that q € Ci\p], and let Uj be a 
direction such that (uj, x) < for all x E Cj. We say that q is j-extremal for p if q = argmax{(p', Uj) \ 
p' G Ci[p] PI P \ {p}}. That is, q is the nearest point to p in this cone, in the (— -direction. Clearly, a 
point p has at most s j-extremal points, one for every admissible cone Ci[p], for any fixed j. See Figure [TOl 
(left). 

(ii) For < i < k, let C[ denote the extended cone that is the union of the seven consecutive cones 
Ci-3, . . . , Cj+3. Let p,q G P, let i be the index such that q G Ci[p], and let Uj be a direction such that 
{uj,x) < for all x G C[ (such Uj's exist if a is smaller than some appropriate constant). We say that the 
point q G P is strongly j-extremal for p if g = arg max{(p', Uj) \ p' G C- [p] PI P \ {p}}. 

(hi) We say that a pair (p, q) £ P x P is (strongly) (j, tj-extremal, for some < j, I < k — 1, if p is 
(strongly) ^-extremal for q and ^ is (strongly) j-extremal for p. 
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Figure 11: Illustration of the setup in Lemma R~Tt the edge e pq is /3-long around v\p], and the "tip" Aa+qa 
of the cone C[q] is empty. 



Lemma 4.1. Let p,q £ P, and let v be a direction such that the edge e pq appears in VD(P) and is (3-long 
around the ray v\p\. Let C[q] be the fi-cone around the ray from q through p. Then (p, v) > {p' , v) for all 



Proof. Refer to Figure [TT] Without loss of generality, we assume that v is the (+x) -direction and that q 
lies above and to the right of p. (In this case the slope of the bisector b pq is negative. Note that q has to 
lie to the right of p, for otherwise b pq would not cross v[p\.) Let v + (resp., v~) be the direction that makes 
a counterclockwise (resp., clockwise) angle of /3 with v. Let a + (resp., a~) be the intersection of e pq with 
v + [p] (resp., with u~[p]); by assumption, both points exist. Let h be the vertical line passing through p. Let 
cr + (resp., a~) be the intersection point of h with the ray emanating from a + (resp., a~) in the direction 
opposite to v~ (resp., v + ); see Figure [TTI 

Note that Apa + a + = 2/3, and that ||a + <7 + || = ||pa + || = ||</a + ||, i.e., a + is the circumcenter of Apa + q. 
Therefore Za + qp = ^Za + a + p = (3. That is, a + is the intersection of the upper ray of C[q] with h. 
Similarly, a~ is the intersection of the lower ray of C[q] with h. Moreover, if there exists a point x € P 
properly inside the triangle Apqa + then ||a + x|| < ||a + p||, contradicting the fact that a + is on e pq . So 
the interior of Apqa + (including the relative interiors of edges pq, a + q) is disjoint from P. Similarly, by 
a symmetric argument, no points of P lie inside Apqa~ or on the relative interiors of its edges pq, a~q. 
Hence, the portion of C[q] to the right of p is openly disjoint from P, and therefore p is a rightmost point of 
P (extreme in the v direction) inside C[q\. □ 

Corollary 4.2. Let p,q € P. (i) If the edge e pq is 3a-long in VD(P) then there are < j,£ < k for which 
(p, q) is a (j, I) -extremal pair. ( ii) If the edge e pq is 9a-long in VD(P) then there are < j, £ < kfor which 
(p, q) is a strongly (j, l)-extremal pair. 

Proof. To prove part (i), choose < j,£ < k, such that e pq is a-long around each of ue\p] and Uj[q], By 
Lemma |4~T1 p is u^-extremal in the a-cone C[q] around the ray from q through p. Let i be the index such 
that p G Ci[q\. Since the opening angle of C[q] is 2a, it follows that Ci[q] C C[q], so p is ^-extremal with 
respect to q, and, symmetrically, q is j-extremal with respect to p. To prove part (ii) choose < j, £ < k, 
such that e pq is 4a-long around each of ug [p] and Uj [q] and apply Lemma |4~T1 as in the proof of part (i). □ 



p'ePn C[q] \ {q}. 
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The stable Delaunay graph. We kinetically maintain a (10a, a)-stable Delaunay graph, whose precise 
definition is given below, using a data-structure which is based on a collection of 2-dimensional orthogonal 
range trees similar to the ones used in flT|. 

Fix < i < s, and choose a "sheared" coordinate frame in which the rays ui and Ui + \ form the x- 
and y-axes, respectively. That is, in this coordinate frame, q £ Cj[p] if and only if q lies in the upper-right 
quadrant anchored at p. 

We define a 2-dimensional range tree Tj consisting of a primary balanced binary search tree with the 
points of P stored at its leaves ordered by their x-coordinates, and of secondary trees, introduced below. 
Each internal node v of the primary tree of Tj is associated with the canonical subset P v of all points that 
are stored at the leaves of the subtree rooted at v. A point p £ P v is said to be red (resp., blue) in P v 
if it is stored at the subtree rooted at the left (resp., right) child of v in Tj. For each primary node v we 
maintain a secondary balanced binary search tree 7j', whose leaves store the points of P v ordered by their 
y-coordinates. We refer to a node w in a secondary tree 7? as a secondary node w o/Tj. 

Each node uiofa secondary tree 7f is associated with a canonical subset P w C P v of points stored at 
the leaves of the subtree of 7% rooted at w. We also associate with w the sets R w C P w and B w C P w of 
points residing in the left (resp., right) subtree of w and are red (resp., blue) in P v . It is easy to verify that 
the sum of the sizes of the sets R w and B w over all secondaiy nodes of Tj is 0(n log 2 n). 

For each secondary node w £ Tj and each < j < k we maintain the points 

= ar S max(p, uj), £?j(w) = arg max(p, uj), 

provided that both R W ,B W are not empty. See Figure [12] It is straightforward to show that if (p, q) is a 
(j, £)-extremal pair, so that q £ C{\p], then there is a secondaiy node w € Tj for which q = Cfj( w ) an d 
p = tf e (w). 




Figure 12: The points £^(i«),£§(tu) for a secondary node w of 

For each p £ P we construct a set containing all points q £ P for which (p, q) is a (J, £)-extremal 
pair, for some pair of indices < j, I < k. Specifically, for each < i < s, and each secondaiy node 
w £ Tj such that p = £fg(w) for some < £ < k, we include in N[p] all the points q such that q = ifj(w) 
for some < j < k. Similarly, for each < i < s, and each secondaiy node w £ Tj such that p = if^{w) 
for some < I < k we include in all the points g such that g = ^{w) for some < j < k. It is easy 
to verify that, for each (i, ^)-extremal pair (p, q), for some < j, £ < k, q is placed in N[p] by the preceding 
process. The converse, however, does not always hold, so in general {p} x N[p] is a superset of the pairs 
that we want. 
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For each < i < s, each point p € P belongs to 0(log 2 n) sets R w and B w , so the size of N[p] is 
bounded by 0(s 2 log 2 n). Indeed, p may be coupled with up to k = 2s neighbors at each of the 0(s log 2 n) 
nodes containing it. 

For each point p € P and < I < k we maintain all points in N[p] in a kinetic and dynamic tournament 
T>i\p\ whose winner q minimizes the directional distance <pe\p, q], as given in (fl]). That is, the winner in 
T>i[p] is Ni[p] in the Voronoi diagram of {p} U N[p]. 

We are now ready to define the stable Delaunay graph G that we maintain. For each pair of points 
p, q G P we add the edge (p, q) to G if the following hold. 

(Gl) There is an index < £ < k such that q wins the 8 consecutive tournaments £>^[p], . . . , T>£ +7 \p\. 

(G2) The point p is strongly {£ + 3) -extremal and strongly (£ + 4) -extremal for q. 

The (10a, a)-stability of G is implied by a combination of Theorems 14. 3 1 and 1 4.41 

Theorem 4.3. For every Wa-long edge e pq € VD(P), the graph G contains the edge (p, q). 

Proof. By Corollary I4.2I (T). there are j and £ such that (p, q) is a (j, £)-extremal pair. By the preceding dis- 
cussion this implies that q is in N[p] . Now since e pq is lOa-long there is an £' such that Np [p] , . . . , Np +7 [p] = 
q in VD(P), and therefore also in the Voronoi diagram of {p} U N[p]. So it follows that q indeed wins the 
tournaments T>(i [p], . . . , V^ +7 [p}. 

By the proof of Corollary 14.21 (h). p is strongly {£' + 3)-extremal and strongly {£' + 4)-extremal for 
q. □ 

Theorem 4.4. For every edge (p, q) € G, the edge e pq belongs to VD(P) and is a-long there. 

Proof. Since (p, q) € G we know that q is in !N[p] and wins the tournaments V^[p], T>£ + \\p], . . . ,T>g +7 [p], 
for some < £ < k and that the point p is strongly (£ + 3)-extremal and {£ + 4)-extremal for q. We prove 
that the rays U£ +3 [p] and n^ +4 [p] stab e pq , from which the theorem follows. 

Assume then that one of the rays U£ + z [p] , ii£+4 \p] does not stab e pq ; suppose it is the ray U£ + 4[q]. (This 
includes the case when e pq is not present at all in VD(P).) By definition, this means that r = A^ +4 [p] ^ q. 
We use Lemma 1431 given shortly below, to show that q cannot win in at least one of the tournaments among 
T>£ [p] , . . . , T>e +7 \p] and thereby get a contradiction. 

According to Lemma 1431 there exists a point r such that ipg+^lp, r] < ip£ + 4\p, q] and p is (£ + 4)- 
extremal for r. Let x = t^+4[p] H b pr and let D be the circle which is centered at x, and passes through r 
and p; see Figure [T3l 

We consider the case where r is to the right of the line from p to q; the other case is treated symmetrically. 
In this case the intersection of b pr and b pq is to the left of the directed line from p to x. Let < i < k — lbe 
the index for which r 6 Cj[p]. If i < s — 1 then there is a secondary node w in the tree Tj for which p £ 
and r € and since p is (£ + 4)-extremal for r, £^ , 4 (w) is equal to p. If i > s then, symmetrically, we 
have a node £ Tj_ s such that r 6 i? w and p E and £,f i+4 (w) is equal to p. We assume that i < s — 1 
in the sequel; the other case is treated in a fully symmetric manner. 

Let v[r] be the ray from r through x, for an appropriate direction v € S 1 , and let Uj be the direction 
which lies counterclockwise to v and forms with it an angle of at least a and at most 2a. Put r' = £,fj{w), 
implying that r' € C;[p] and (r',Uj) > (r,Uj). In particular, r' belongs to N[p]. If r' is inside D (and in 
particular if r 1 = r) then q cannot win the tournament ^+4[p] which is the contradiction we are after. So 
we may assume that r' is outside D. 
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Figure 13: Proof of Theorem l4~4l the case when r is to the right of the line from p to q. The line h orthogonal 
to Uj through r intersects the circle D at a point y outside Ci[p], which implies that r' is to the right of the 
line from p to r. Assuming r ^ r', the point z = b pr n b pr < is inside the cone bounded by u i+i [p] and ue +7 \p\. 
Hence, u e+7 [p] hits b pr > before b pr . 

Let h be the line through r orthogonal to Uj-. Clearly, h intersects D at two points, r and another point 
y (lying counterclockwise to r along 3D, by the choice of Uj). Since Zrpy = \/-Txy, and Zrxy equal to 
twice the angle between v and Uj, Zrpy is at least a, so y is outside Ci[p\. By assumption, r' lies in the 
halfplane bounded by h and containing p. Since we assume that r' is not in D it must be to the right of the 
line from x to r. It follows that b pr i intersects b pr at some point z to the right of the line from p to x; see 
Figure [13] 

We claim that z is inside the cone with apex p bounded by the rays [p] and U£ + j [p] . Indeed, suppose 
to the contrary that the claim is false. It follows that in the diagram VD({r, r',p}) the edge e pr is a- 
long around Uj[r]. Indeed, denote the intersection point of U£ +7 [p] and b pr as w (see Figure [T3"T). Then 
Zxrw = Zxpw = 3a. Since the angle between v[r] and Uj[r] is between a and 2a, the claim follows. 
Now, according to Lemma |4~T1 (r,Uj) < {r',Uj}, which contradicts the choice of r'. It follows that z is 
in the cone bounded by ue+^lp] and u^ + j[p] and thus ue + j[p] hits b pr > before 6 pr , and therefore also before 
b pq . Hence, q cannot win T>e +7 \p], and we get the final contradiction which completes the proof of the 
theorem. □ 

Remark: We have not made any serious attempt to reduce the constants c appearing in the definitions of 
various (ca, a)-SDGs that we maintain. We suspect, though, that they can be significantly reduced. 
To complete the proof of Theorem 14.41 we provide the missing lemma. 

Lemma 4.5. Let p,q € P be a pair of points and < £ < k — 1 an index, such that the point p is strongly 
t- extremal for q but Np[p] ^ q. Then there exists a point r such that <pi\p, r] < fe\p, q] and p is (.-extremal 
for r. 

Proof. Let < i < k — lbe the index for which q 6 C-i [p] and let h be the line through p, orthogonal to U£. 
Assume without loss of generality that h is vertical and the ray ui [p] extends to the right of h. 

Let a be the point at which U(\p] intersects the bisector b pq , and let D be the disk centered at a whose 
boundary contains both p and q. Since Ng\p] ^ q, the interior of D must contain some other point r € P; 
see Figure [14] 
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Let C[q] be the cone emanating from q such that each of its bounding rays makes an angle of (3 = 3a 
with the ray from q through p; in particular C[q] contains p. Let a + (resp., a~) denote the upper (resp., 
lower) endpoint of the intersection of C[q) and h. Since p is strongly ^-extremal for q, the interior of the 
triangle Aa + qa~ does not contain any points of P. Hence, r must be outside the triangle Aa + qa~. So 
either r is above qa + (and inside D) or below qo~ (and inside D). 

Assume, without loss of generality, that r is below qa~ , as shown in Figure [14] (The case where r is 
above qa + is fully symmetric.) Let t + and t~ denote the intersection points qa + n 3D and qo~ n 3D, 
respectively. Let e be the point at which the ray from r through t~ intersects h. Then the intersection of the 
triangle Aa + re and Aa + qa~ is empty. Among all the points of P in D we choose r so that its x-coordinate 
is the smallest. For this choice of r we also have that Aa+re \ Ao + qa~ is empty (since it is contained in 
D and lies to the left of r). In other words, Aa + re is empty. 

Let 7 + (resp., 7"") denote the angle Zpra + (resp., Z.prt~). It remains to show that 7 + > |/3 and 
7~ > |/3. This will imply that the cone CV [r] that contains p is fully contained in the cone bounded by the 
rays from r through a + and t~, so p is extreme in the u^-direction within Cy [r], which is what the lemma 
asserts. Since r is inside D, it is clear that 7" > Zpqt~ = f3. The angle 7 + however may be smaller than 
ft, but, as we next show, tan 7+ > | tan j3. Indeed, fix an angle 6 and let T(6) denote the circular arc which 
is the locus of all points z that are to the right of h and the angle Zpza + is 9. The endpoints of T(0) are p 
and cr + , and its center a* is on the (horizontal) bisector of pa + ; see Figure [T51 (left). 

Notice that T{6) intersects dD at two points, one of which is p, which are symmetric with respect to the 
line through a and a*. As 6 decreases a* moves to the right, and the intersection of T(9) with dD rotates 
clockwise around dD. Consider the smallest 9 such that T(0) intersects D on or below qt~. It follows that 
this intersection is at t~ . See Figure [T51 (right). 

This shows that for fixed p and q, the position of r in D below the line qt~ which minimizes 7+ is at 
t~. To complete the analysis, we look for the position of q that minimizes 7+ when r is at t~. Note that, as 
q moves along dD, the points t + and t~~ do not change. As shown in Figure [T6l (left). 7 + decreases when 
q tends counterclockwise to t + . When q is at t + , qa + is tangent to D. A simple calculation, illustrated 
in Figure [TBI (right), shows that tan 7+ = |tan/3. By the inequality tan(3x) > 3tanx, for x sufficiently 
small, it follows that 7+ > i/3, implying, as noted above, that the point p is ^-extremal for r. This completes 
the proof of the lemma. □ 
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Figure 15: Left: The circular arc T(6) is the locus of all points which are to the right of p<r + and see it at angle 
0. Right: To minimize 8 we increase the radius of T{9) until one of its intersection points with D coincides 
with t~~ . 




In Section |4~T1 we describe a naive algorithm for kinetic maintenance of G, which encounters a total of 
0*(k A n 2 ) events in the tournaments T>g\p]. In Section 1421 we consider a slightly more economical definition 
of the tournaments T>i\p\, yielding a solution which processes only 0*(k 2 n 2 ) events in 0*(k 2 n 2 ) overall 
time. 

4.1 Naive maintenance of G 

As the points of P move, we need to update the SDG G, which, as we recall, contains those edges (p, q) 
such that q wins 8 consecutive tournaments Pf[p], ... ,T)£ + j[p] of p, and p is strongly (£ + 3)-extremal 
and (I + 4) -extremal for q. We thus need to detect and process instances at which one of these conditions 
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changes. There are several events at which such a change can occur: 

(a) A change in the sets of neighbors N[p] , for p ^ P. 

(b) A change in the status of being strongly ^-extremal for some pair (p, q). 

(c) A change in the winner of some tournament Vg[p] (at which two existing members of attain the 
same minimum distance in the direction ug). 

Note that each of the events (a)-(b) can arise only during a swap of two points in one of the s directions 
Uq, . . . , u s -i or in one of the directions orthogonal to these vectors. 

For each < i < s - 1 we maintain two lists. The first list, Lj, stores the points of P ordered by their 
projections on a line in the Uj-direction, and the second list, K t , stores the points ordered by their projections 
on a line orthogonal to the m -direction. We note that, as long as the order in each of the 2s lists Ki, Li 
remains unchanged, the discrete structure of the range trees Tj, and the auxiliary items ^fi{vj),^fj(w), does 
not change either. More precisely, the structure of % changes only when two consecutive elements in Ki 
or in Ki + \ swap their order in the respective list; whereas the auxiliary items £,^ t j(w),£ : fj(w), stored at 
secondary nodes of Tj, may also change when two consecutive points swap their order in the list Lj. There 
are 0(sn 2 ) = 0(n 2 ) discrete events where consecutive points in Ki or Lj swap. We call these events K{- 
swaps and Li-swaps, respectively. Each such event happens when the line trough a pair of points becomes 
orthogonal or parallel to Uj. We can maintain each list in linear space for a total of O(sn) space for all lists. 
Processing a swap takes O(logn) time to replace a constant number of elements in the event queue (and 
more time to update the various structures, as discussed next). 

The range trees Tj. As just noted, the structure of Tj changes either at a fQ-swap or at a .fQ+i-swap. As 
described in [fl] Section 4], we can update Tj when such a swap occurs, including the various auxiliary data 
that it stores, in 0(s log 2 n) time. (The factor s is due to the fact that we maintain O(s) extreme points 
£,fl{w) and if"j{w) in each secondary node w of Tj, whereas in [1] only two points are maintained.) 

In a similar manner, an Lj-swap of two points p, q may affect one of the items £^ (w) and (w) stored 
at any secondary node w of any Tj, for < i < s — 1, such that both p, q belong to R w or to B w . Each 
Tj has only 0(log 2 n) such nodes, and the data structure of [fl] allows us to update Tj, when an Lj-swap 
occurs in 0(log 2 n) time. Summing up over all < i < s — 1, we get that the total update time of the 
range trees after an Lj-swap is 0(s log 2 n). As follows from the analysis in UJ Section 4], the trees Tj, for 
< i < s — 1, require a total of 0{s 2 n log n) storage (because of the O(s) items ifi{w),if'j{w) stored at 
each secondary node of each of the s trees). 

The tournaments T)g\p]. The kinetic tournament T>c[p], for p € P and < I < k — 1 contains the points in 
the set Since 3sf[p] varies both kinetically and dynamically and therefore the tournaments Vi[p] need 
to be maintained as kinetic and dynamic tournaments, in the manner reviewed in Section [2] 

For < i < s — 1, we define II to be the set of pairs of points (p, q), such that there exists a secondary 
node iu in 7i, and indices < j,£ < k — 1, for which p = £fg(w) and q = £fj(w). For a fixed i, a point 
p belongs to 0(s log 2 n) pairs (p, q) in Ilj, for a total of 0(s 2 log 2 n) pairs over all sets Ilj. It follows that 
the total size of all the sets Ilj is 0(s 2 n log 2 n). Any secondary node of any tree Tj, for < i < s — 1, 
contributes at most 0(s 2 ) pairs to the respective set Ilj. 

The set N[p] consists of all the points q such that there exists a set Ilj that contains the pair (p, q) or the 
pair (q,p). So the total size of the sets over all points p, is 0(s 2 n log 2 n). A set N[p] changes only 
when one of the sets Ilj changes, which can happen only as the result of a swap. 

Specifically, when ^(w) changes for some < i < s — 1 and < £ < k — 1, from a point p to a point 
p', we make the following updates, (i) If p ^ (w) for all £' ^ £ then for every < j < k — 1 we delete 
the pair (p, A(w)) from Ilj. (ii) We add the pair (p f , ifj{w)) to Ilj. We make analogous updates when one 
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of the values £fj{w) changes. When a node w is created, deleted, or involved in a rotation, we update the 
pairs iifi(w), £f"Aw)) in Ilj for every £ and j. In such a case we say that node w is changed. 

A change of ^(w) or £fj(w) in an existing node w generates O(s) changes in Ilj and thereby O(s) 
changes to the sets N[p]. Thus, it may generate 0(s 2 ) updates to the tournaments Z^[p]. A change of a 
secondary node may generate 0(s 2 ) changes to the sets 3Sf[p] and thereby 0(s 3 ) updates to the tournaments 

A point £fg(w) or £fi(w) changes during either a Ki, Ki+\, or L^-swap. Each L^-swap, for any £, 
causes 0(slog 2 n) points ^(w) or £fi(w) to change (over the entire collection of trees), and therefore 
each swap causes 0(s 3 log 2 n) updates to the tournamnets T>#\p]. The number of nodes which change in 
Tj by a Ki or ifj + i-swap is 0(log 2 n). Each such change causes 0(s 3 ) updates to the tournaments T>i[p\. 
Therefore the total number of updates to tournaments due to changes of nodes is also 0(s 3 log 2 n) per swap. 

The number of swaps is 0(sn 2 ), so overall we get 0{s A n 2 log 2 n) updates to the tournaments. The size 
of each individual tournament is 0(s 2 log 2 n). By Theorem 12.11 these updates generate 

0(s 4 n 2 log 2 n ■ /3 r+ 2(s 2 log 2 n) log(s 2 log 2 n)) = 0(s 4 n 2 /3 r +2(s log n) log 2 n log(s log n)) 

tournament events, which are processed in 

0(s i n 2 log 2 n ■ f3 r+ 2(s 2 log 2 n) log 2 (s 2 log 2 n)) = 0(s 4 n 2 ■ /3 r +2(s log n) log 2 n log 2 (s log n)) 

time. Processing each individual tournament event takes (9(log 2 log n + log 2 s) time. 

Since the size of each tournament is 0(s 2 log 2 n) and there are 0(ns) tournaments, the total size of all 
tournaments is 0(s 3 n log 2 n). 

Testing whether p is strongly ^-extremal for the winner of 2^[p]. For each < i < s — 1, and for each 
pair (p, q) 6 Ilj we maintain those indices < t < k — 1 (if there are any) for which p is strongly ^-extremal 
for q. Recall that each point p belongs to 0(s 2 log 2 n) pairs in the sets Ilj. 

We use the trees Jj for i — 3 < j < i + 3 to find, for a query q, the point arg m.ax q > e p n c'.[q] (<?'> u t)> f° r 
each < I < k — 1. The query time is 0(s log 2 n) Using this information we easily determine, for a pair 
(p, q), for which values of I p is strongly ^-extremal for q. 

As explained above, every swap changes 0(s 2 log 2 n) pairs of the sets Ilj. When a new pair is added to 
a set Ilj we query the trees Tj, i — 3 < j < i + 3, to find for which values of £, p is strongly ^-extremal for 
q (and vice versa). This takes a total of 0(s 3 log 4 n) time for each swap. 

Furthermore, a point p can cease (or start) being strongly ^-extremal for q only during a swap which 
involves either p or q. So when we process a swap between p and some other point we recompute, for all 
pairs (p, x) and (x,p) in the current sets Ilj and for every < £ < k — 1, whether p is strongly ^-extremal 
for x, and whether x remains strongly ^-extremal for p. This adds an overhead of 0(s 3 log 4 n) time at each 
swap. 

The following theorem summarizes the results obtained so far in this section. 

Theorem 4.6. The SDG G can be maintained using a data structure which requires O ((n/a 3 ) log 2 n) 
space and encounters two types of events: swaps and tournament events. 
There are 0(n 2 /a) swaps, each processed in O (log 4 n/a 3 ) time. There are 

O ((n 2 /a 4 ) log 2 n/3 r+ 2 (log n/a) log (log n/a)) 
tournament events which are processed in overall 

O ((n 2 /a 4 ) log 2 n(3 r+2 (log n/a) log 2 (log n/a)) 
time. Processing each individual tournament event takes 0(log 2 log n + log 2 (l/a)) time. 
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4.2 An even faster data structure 

We next reduce the overall time and space required to maintain G roughly by factors of s 2 and s, respectively 
(bringing the dependence on s of both bounds down to roughly s 2 ). We achieve that by restricting each 
tournament T>g\p] to contain a carefully chosen subset C of size 0(s log 2 n) (recall that the size 
of the entire set N[p] is 0(s 2 log 2 n)). The definition of 3sf^[p] is based on the following lemma. Its simple 
proof is given in Figure ITTI 

Lemma 4.7. Let p,q € P and let i be the index for which q £ Ci\p\. Let < £ < k — 1 be an index, and 
v € S 1 a direction such that the rays U(\p\ and v[q] intersect b pg at the same point. Then v lies in one of the 
two consecutive cones C^^, C^u(\ +1 , where ((i,£) = 2i + s — t. 




Figure 17: Proof of Lemma 14771 We assume that q e Cf[p], and that the rays u e [p] and v[q] hit b pq at the 
same point w. Then the angle x = Zwpq = (i + 1 - £)a - t, for some < t < a. The orientation of qp is 
(i + l)a - 1 + 7r = (i + s + 1)q - 1. Hence, the orientation of v is (i + s + l)a — t + x= (2i + s — £ + 2)a - 2t. 
Thus, the direction v lies in the union of the two consecutive cones C c(i ^), C c(i ^ +1 , for C(M) = 2i + s - £. 

It follows that in Corollary 14.21 we can require that the indices < j, I < k — 1, for which (p, q) is a 
(strongly) (j, -extremal pair, satisfy ((i,£) < j < (,{i,£) + 2. Indeed, we may require that the vectors 
Uj -,[q], U([p] hit b pq at the respective points x and y for which the angle Zxpy = Zxqy is at most a, which, 
in turn, happens only if Uj bounds one of the cones C^an, C^u(\ + x. 

For all < i < s — 1 and < £ < k — lwe define a set IL^ which consists of all pairs (p, q) of 
points of P such that there exists a secondary node w in Tj, and indices £ and C(h^) < 3 < C(h^) + 2, 
such that p = ^(w) and (/ = Cfj( w ) or p = ^(w) and g = C^jiw). We define the set N^[p] to consist of 
all points g such that (p, g) € Ilj ^. For a point p the set of points that participate in the reduced tournament 
V(\p\ is \p\. (Note that this rule distributes a point q G Ng[p] to only seven nearby tournaments. 

Nevertheless, when the edge pq is sufficiently long, q will belong to several consecutive neighborhoods 
and therefore will appeal - in more tournaments, in particular in at least eight consecutive tournaments 
at which it should win, according to the definition of our SDG.) 

We claim that, with this redefinition of the tournaments T>(\p], Theorems I4.3l and [4.4l still hold. To verify 
that Theorem 14. 3 I holds one has to follow its (short) proof and notice that, by Lemma l4~7l the point q belongs 
to the eight reduced tournaments which it is supposed to win. 

We next indicate the changes required in the proof of Theorem 14.41 We use the same notation as in the 
original proof of Theorem 14.41 and recall that it assumed by contradiction that, say, Ni + i\p] ^ q even though 
q wins the tournaments T>i\p],T>£ +1 [p], . . . , D^ +7 [p], and the point p is strongly (^+3)- and (£+4) -extremal 
for q. We use Lemma 1431 to establish the existence of some point r 6 P such that 9^+4 [p, r] < </^ + 4[p, q] 
and p is (£ + 4) -extremal for r. Let i be the index for which r € C-i [p], and let w be the secondary node in 
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Tj for which r G B w and p G R w . Note that p = £f^ +4 (w). We next cnoose an index j such that the point 
r' = C^jiw) either satisfies that ipi +7 {p, r'] < tp£ +7 \p, q] if r is to the right of the line from p to q, or that 
(fe+i [p, r'] < (fi+i \p, q] if r is to the left of the line from p to q. To re-establish Theorem 14.4 l it suffices to 
show that r' participates in the reduced tournament Vi +7 [p} (resp., if r is to the right (resp., left) 

of the line from p to q. 

It follows from the way we denned j in the original proof and from Lemma 14771 that £ + 4) — 2 < 
j < ((i,£ + 4:)-l (if r is to the right of the line from p to g) or ((i,£ + 4) + l <j< ((i, £ + 4) + 2 (if r is to 
the left of the line from p to q). So r' G N£+4[p] and therefore r' does participate in the reduced tournament 
T>£ + i [p] or T>£ +7 [p] . Indeed, the direction v used in that proof lies in one of the cones C^u^ + ^\ , CW^+^+i . 
The direction uj then fomis an angle between a and 2a with v, which lies counterclockwise from v if r lies 
to the right of the line from p to q, or clockwise from v in the other case. This is easily seen to imply the 
two corresponding constraints on j; see Figure [T3l 

We change our algorithm accordingly to maintain only the reduced tournaments. 

Now every secondary node w of any range tree Tj contributes only seven pairs to each set IL^, for 
< £ < k — 1, so the size of each such set is O(nlogn). Since there are 0(s 2 ) sets Hij, their total size 
is 0(s 2 n log n). Each pair in each Ilj^ contributes an item to a constant number of tournaments, so the 
total size of the tournaments is 0(s 2 nlogn). Each individual tournament V^[p\ is now of size 0(s log 2 n), 
because p belongs to 0(log 2 n) pairs in each set 11^/ for < i < s — 1, < £' < k — 1, and T>g\p\ inherits 
only those points q that come from pairs (p, q) G Hi p, f or < i < s — 1 and £ — 3 <£'<£ + 3. 

When £^ [tu] changes from p to p' for some < i < s — 1 and < £ < k — 1, at most a constant number 
of pairs (p,^j(w)) for ((i,£) < j < C{i,£) + 2 are deleted from H^, and a constant number of pairs 
(p'j (' u; )) f° r C(^) ^) < J < ^) + 2 are added to Ilj^. Similar changes take place in Hij for those three 
indices j satisfying < £ < C(hj) + 2- When C/jfH changes from g to q' for some < i < s — 1 

and < j < k — 1, at most a constant number of pairs {if^{w), q) are deleted from Iljj for the indices 
^ satisfying C(hj) < ^ < C(hj) + 2, and a constant number of pairs {£,fn{w),q') are added for the same 
values of £. Similarly, at most a constant number of pairs {£,fi{w),q) are deleted from H ^ for the indices 
£ satisfying ((i,£) < j < ((i,£) + 2, and a constant number of pairs (£fi(w), q') are added for the same 
values of £. 

A change of a secondary node w in the tree causes O(s) pairs in the sets IL^ to change. 

Any Kj-swap changes 0(log 2 n) nodes in Tj and thereby causes 0(s log 2 n) pairs in the sets IL;^ to 
change. Any Lj-swap changes 0(s log 2 n) extremal points ^f - [w], ^ [w] at secondary nodes w of the trees 
Tj, and thereby causes 0(slog 2 n) pairs in the sets IL;^ to change. Since each pair in TL^ contributes an 
item to a constant number of tournaments it follows that 0(s log 2 n) points are inserted to and deleted from 
the tournaments Ug\p] at each swap. 

According to Theorem 12. li the size of each tournament is 0(s log 2 n) - the number of elements that it 
contains. So the total size of all tournaments is (9(s 2 nlogn). In total we get that there are 0{s 2 n 2 log 2 n) 
updates to tournaments during swaps. These updates generate 

0(s 2 n 2 log 2 n/3 r+ 2(s logn) log(s log n)) 

tournament events that are processed in overall 

0(s 2 n 2 log 2 nf3 r+ 2(s log n) log 2 (s logn)) 

time. Each individual tournament event is processed in 0(log 2 logn + log 2 s) time and each swap can be 
processed in 0(s log 2 n log 2 (s log n)) time. 
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In addition, for each pair (p, q) € 11^ we record whether p is strongly ^-extremal for q. We maintain 
this information using the trees Tj, for i — 3 < j < i + 3, as described above, which allow for any p,q € P 
and < I < k — lto test, in 0(log 2 n) time, if p is strongly ^-extremal for q. At each swap event we spend 
0(s log 4 n) extra time to compute for 0(s log 2 n) pairs (p, q) which are added to the sets whether p is 
strongly ^-extremal for q. 

Consider a pair (p, q) € JJj^. The point p may stop being strongly ^-extremal for q only during a swap 
which involves p or q. So, as before, at each swap we find the 0(s log 2 n) pairs containing one of the points 
involved in the swap, and recompute, in 0(s log 4 n) total time, for each such pair (p, q), whether the strong 
extremal relation holds. We thus obtain the following summary result. 

Theorem 4.8. Let P be a set of n moving points in M 2 under algebraic motion of bounded degree, and let 
a > Obe a sufficiently small parameter. A (10a, a)-SDG ofP can be maintained using a data structure that 
requires 0((n/a 2 ) logn) space and encounters two types of events: swap events and tournament events. 
There are 0(n 2 /a) swap events, each processed in 0(log 4 (n)/a) time. There are 

0{(n/a) 2 f3 r+2 {log{n)/a) log 2 nlog(log(n)/a)) 

tournament events, which are handled in a total of 

O ( (n/a) 2 f3 r+2 (log (n) /a) log 2 n log 2 (log (n) /a) ) 

processing time. The worst-case processing time of a tournament event is 0(log 2 (log(n)/a)). The data 
structure is also local, in the sense that each point is stored, at any given time, at only 0(log 2 n/a 2 ) places 
in the structure. 

Concerning locality, we note that a point participates in O(s) projection tournaments at each of 0(s log 2 n) 
tree nodes. If it wins in at least one of the projection tournaments at a node, it is fed to O(s) directional 
tournaments. So it appears in 0(s 2 log n) places. 

Remarks: (1) Comparing this algorithm with the space-inefficient one of Section 12 we note that they both 
use the same kind of tournaments, but here much fewer pairs of points (0*(n/a 2 ) instead of 0(n 2 /a)) 
participate in the tournaments. The price we have to pay is that the test for an edge pq to be stable is 
more involved. Moreover, keeping track of the subset of pairs that participate in the tournaments requires 
additional work, which is facilitated by the range trees Tj. 

(2) To be fair, we note that our 0*(-) notation hides polylogarithmic factors in n. Hence, comparing the anal- 
ysis in this section with Theorem l3.6l we gain when n is smaller than some threshold, which is exponential 
in 1/a. 

5 Properties of SDG 

We conclude the paper by establishing some of the properties of stable Delaunay graphs. 

Near cocircularities do not show up in an SDG. Consider a critical event during the kinetic maintenance 
of the full Delaunay triangulation, in which four points o, b, c, d become cocircular, in this order, along their 
circumcircle, with this circle being empty. Just before the critical event, the Delaunay triangulation involved 
two triangles, say, abc, acd. The Voronoi edge e ac shrinks to a point (namely, to the circumcenter of abed 
at the critical event), and, after the critical cocircularity, is replaced by the Voronoi edge eu, which expands 
from the circumcenter as time progresses. 
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Figure 18: The near collinearity that corresponds to a Voronoi edge becoming a-short. 



Our algorithm will detect the possibility of such an event before the criticality occurs, when e ac becomes 
a-short (or even before this happens). It will then remove this edge from the stable subgraph, so the actual 
cocircularity will not be recorded. The new edge will then be detected by the algorithm only when it 
becomes sufficiently long (if this happens at all), and will then enter the stable Delaunay graph. In short, 
critical cocircularities do not arise at all in our scheme. 

As noted in the introduction, a Delaunay edge ab (interior to the hull) is just about to become a-short 
or a-long when the sum of the opposite angles in its two adjacent Delaunay triangles is it — a (see Figure 
U}. This shows that changes in the stable Delaunay graph occur when the "cocircularity defect" of a nearly 
cocircular quadruple (i.e., the difference between ir and the sum of opposite angles in the quadrilateral 
spanned by the quadruple) is between a and ca, where c is the constant used in our definitions in Section 
[3] or Section [4] Note that a degenerate case of cocircularity is a collinearity on the convex hull. Such 
collinearities also do not show up in the stable Delaunay graphJl A hull collinearity between three nodes 
a, b, c is detected before it happens, when (or before) the corresponding Voronoi edge becomes a-short, in 
which case the angle Zacb, where c is the middle point of the (near-)collinearity becomes ir — a (see Figure 
[T~8l >. Therefore a hull edge is removed from the SDG if the Delaunay triangle is almost collinear. The edge 
(or any new edge about to replace it) re-appears in the SDG when its corresponding Voronoi edge is long 
enough, as before. 

SDGs are not too sparse. Consider the Voronoi cell Vor(p) of a point p, and suppose that p has only one 
a-long edge e pg . Since the angle at which p sees e pq is at most ir, the sum of the angles at which p sees the 
other edges is at least n, so Vor(p) has at least ir /a a-short edges. Let m\ denote the number of points p 
with this property. Then the sum of their degrees in DT(P) is at least mi(7r/a + 1). Similarly, if mo points 
do not have any a-long Voronoi edge, then the sum of their degrees is at least 27rmo / a. Any other point at 
least two a-long Voronoi edges and its degree is at least 3 if it is an interior point, or at least 2 otherwise. So 
the number of a-long edges is at least (recall that each a-long edge is counted twice) 

n — mi — tjiq + mi/2 = n — (mi + 2mo)/2. (2) 

Let h denote the number of hull vertices. Since the sum of the degrees is 6n — 2h — 6, we get 

/ 7T \ 7T 

3(n — h — mi — mo) + 2h + mi I hi) + 2mo— < 6n — 2h — 6, 

Va / a 

3 Even if they did show up, no real damage would be done, because the number of such collinearities is only O* (n 2 ); see, e.g., 
(22) 
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Figure 19: If the points of P lie on a sufficiently spaced shifted grid then the number of a-long edges in 
VD(P) (the vertical ones) is close to n. 



implying that 




Plugging this inequality in ©, we conclude that the number of a-long edges is at least 



2(vr/a - 2) 



As a decreases, the number of edges in the SDG is always at least a quantity that gets closer to n. This is 
nearly tight, since there exist n-point sets for which the number of stable edges is only roughly n, see Figure 



Closest pairs, crusts, /3-skeleta, and the SDG. Let (3 > 1, and let P be a set of n points in the plane. The 
(3-skeleton of P is a graph on P that consists of all the edges pq such that the union of the two disks of radius 
(/3/2)d(p, q), touching p and q, does not contain any point of P \ {p, q}. See, e.g., 0[T8l for properties of 
the /3-skeleton, and for its applications in surface reconstruction. We show that the edges of the /3-skeleton 
are a- stable in DT(P), provided (3 > 1 + ft(a 2 ). In Figure [20] we sketch a straightforward proof of the fact 
that the edges of the /3-skeleton are a-stable in DT(P), provided that (3 > 1 + Q(a 2 ). 



Figure 20: An edge pq of the /3-skeleton of P (for (3 > 1). c x and c 2 are centers of the two P-empty disks 
of radius (f3/2)d(p,q) touching p and q. Clearly, each of p, q sees the Voronoi edge e pq at an angle at least 
29 = Zc lP q + Zc 2 pq (so it is 26>-stable). We have 1//3 = cos0 « 1 - 6 2 /2 or /3 = 1 + Q(9 2 ). That is, for 
(3 > 1 + Q(a 2 ) every edge of the ,3-skeleton is a-stable. 

A similar argument shows that the stable Delaunay graph contains the closest pair in P{t) as well as the 
crust of a set of points sampled sufficiently densely along a 1 -dimensional curve (see 0IH for the definition 
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Figure 21: ab is an edge of the relative neighborhood graph but not of SDG. 




Figure 22: A wheel-like configuration that disconnects p in the stable Delaunay graph. The Voronoi diagram 
is drawn with dashed lines, the stable Delaunay edges are drawn as solid, and the remaining Delaunay 
edges as dotted edges. The points of the "wheel" need not be cocircular. 



of crusts and their applications in surface reconstruction). We only sketch the argument for closest pairs: If 
(p, q) is a closest pair then pq € DT(P), and the two adjacent Delaunay triangles Apqr + , Apqr~ are such 
that their angles of r + , r~ are at most 7r/3 each, so e pq is (7r/3)-long, ensuring that pq belongs to any stable 
subgraph for a sufficiently small; see Q for more details. We omit the proof for crusts, which is fairly 
straightforward. 

In contrast, stable Delaunay graphs need not contain all the edges of several other important subgraphs of 
the Delaunay triangulation, including the Euclidean minimum spanning tree, the Gabriel graph, the relative 
neighborhood graph, and the all-nearest-neighbors graph. An illustration for the relative neighborhood 
graph is given in Figure [21] As a matter of fact, the stable Delaunay graph need not even be connected, as is 
illustrated in Figure |22] 

Completing SDG into a triangulation. As argued above, the Delaunay edges that are missing in the 
stable subgraph correspond to nearly cocircular quadruples of points, or to nearly collinear triples of points 
near the boundary of the convex hull. Arguably, these missing edges carry little information, because they 
may "flicker" in and out of the Delaunay triangulation even when the points move just slightly (so that all 
angles determined by the triples of points change only slightly). Nevertheless, in many applications it is 
desirable (or essential) to complete the stable subgraph into some triangulation, preferrably one that is also 
stable in the combinatorial sense — it undergoes only nearly quadratically many topological changes. 

By the analysis in Section [3] we can achieve part of this goal by maintaining the full Delaunay triangu- 
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Figure 23: The triangulation DT°(P) of an 8-point set P. The points a, b, c, d, which do not lie on the convex 
hull of P, still lie on the boundary of the union of the triangles of DT°(P) because, for each of these points 
we can place an arbitrary large homothetic interior-empty copy of Q which touches that point. 



lation DT°(P) under the polygonal norm induced by the regular /c-gon Q^. This diagram experiences only 
a nearly quadratic number of topological changes, is easy to maintain, and contains all the stable Euclidean 
Delaunay edges, for an appropriate choice of k ~ 1/a. Moreover, the union of its triangles is simply con- 
nected — it has no holes. Unfortunately, in general it is not a triangulation of the entire convex hull of P, as 
illustrated in Figure [23] 

For the time being, we leave it as an open problem to come up with a simple and "stable" scheme for 
filling the gaps between the triangles of DT°(P) and the edges of the convex hull. It might be possible to 
extend the kinetic triangulation scheme developed in 11161 . so as to kinetically maintain a triangulation of 
the "fringes" between DT°(i- > ) and the convex hull of P, which is simple to define, easy to maintain, and 
undergoes only nearly quadratically many topological changes. 

Of course, if we only want to maintain a triangulation of P that experiences only a nearly quadratically 
many topological changes, then we can use the scheme in lPT6ll . or the earlier, somewhat more involved 
scheme in [2j. However, if we want to keep the triangulation "as Delaunay as possible", we should include 
in it the stable portion of DT, and then the efficient completion of it, as mentioned above, becomes an issue, 
not yet resolved. 

Nearly Euclidean norms and some of their properties. One way of interpreting the results of Section 
3 is that the stability of Delaunay edges is preserved, in an appropriately defined sense, if we replace the 
Euclidean norm by the polygonal norm induced by the regular /c-gon Qfc (f° r k « 1/a). That is, stable 
edges in one Delaunay triangulation are also edges of the other triangulation, and are stable there too. Here 
we note that there is nothing special about The same property holds if we replace the Euclidean norm 
by any sufficiently close norm (or convex distance function |[8l). 

Specifically, let Q be a closed convex set in the plane that is contained in the unit disk Dq and contains 
the disk D' = (cos o)Dq that is concentric with Do and scaled by the factor cos a. This is equivalent to 
requiring that the Hausdorff distance H (Q, Dq) between Q and Do be at most 1 — cos a sa a 2 /2. We define 
the center of Q to coincide with the common center of Dq and D' . 

Q induces a convex distance function cIq, defined by cLq{x : y) = min{A | y G x + XQ}. Consider the 
Voronoi diagram Vor^(P) of P induced by cIq, and the corresponding Delaunay triangulation DT^(P). 
We omit here the detailed analysis of the structure of these diagrams, which is similar to that for the norm 
induced by Qj,, as presented in Section [3] See also EUl for more details. Call an edge e pq of Vor^(P) 
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a-stable if the following property holds: Let u and v be the endpoints of e pq , and let Q u , Q v be the two 
homothetic copies of Q that are centered at u, v, respectively, and touch p and q. Then we require that 
the angle between the supporting lines at p (for simplicity, assume that Q is smooth, and so has a unique 
supporting line at p (and at q); otherwise, the condition should hold for any pair of supporting lines at p or at 
q) to Q u and Q v is at least a, and that the same holds at q. In this case we refer to the edge pq of DT^ (P) 
as a-stable. 

Note that Q^-stability was (implicitly) defined in a different manner in Section [3] based on the number 
of breakpoints of the corresponding Voronoi edges. Nevertheless, it is easy to verify that the two definitions 
are essentially identical. 




Figure 24: An Illustration for Claim 15.1 1 

A useful property of such a set Q is the following: 

Claim 5.1. Let a be a point on dQ and let I be a supporting line to Q at a. Let b be the point on 8Dq closest 
to a (a and b lie on the same radius from the center o). Let 7 be the arc o/ODq, containing b, and bounded 
by the intersection points of £ with ODq. Then the angle between £ and the tangent, r, to Dq at any point 
along 7, is at most a. 

Proof. Denote this angle by 6. Clearly 9 is maximized when r is tangent to Dq at an intersection of £ and 
8Dq. See Figure [24] It is easy to verify that the distance from o to £ is cos 9. But this distance has to be 
at least cos a, or else dQ would have contained a point inside D' , contrary to assumption. Hence we have 
cos 6 > cos a, and thus 9 < a, as claimed. □ 

We need a few more properties: 

Claim 5.2. Let Q\ and Qi be two homothetic copies of Q and let w be a point such that (i) w lies on dQ\ 
and on dQ% and (ii) w and the respective centers o\, 02 of Q\, Q2 are collinear. Then Q\ and Q2 are 
tangent to each other at w; more precisely, they have a common supporting line at w, and, assuming dQ to 
be smooth, w is the only point of intersection of dQ\ PI 8Q2 (otherwise, dQ\ PI 8Q2 is a single connected 
arc containing w.). 

Proof. Map each of Q\, Q2 back to the standard placement of Q, by translation and scaling, and note that 
both transformations map w to the same point wq on dQ. Let £q be a supporting line of Q at wq, and let l\, 
£2 be the forward images of £ under the mappings of Q to Q\ and to Q2, respectively. Clearly, £\ and £2 
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Figure 25: An Illustration for Claim 15.21 



coincide, and are a common supporting line of Q\ and Q2 at w. See Figure [25] The other asserted property 
follows immediately if dQ is smooth, and can easily be shown to hold in the non-smooth case too; we omit 
the routine argument. □ 

Claim 5.3. Let a and b be two points on dQ, and let i a and lb be supporting lines of Q at a and b, 
respectively. Then the difference between the angles that l a and inform with ab is at most 2a. 

Proof. Denote the two angles in the claim by 9 a and #5, respectively. Let a' (resp., b') be the point on ODq 
nearest to (and co-radial with) a (resp., b). Let t\, T2 denote the respective tangents to Dq at a' and at b'. 
Clearly, the respective angles 9\, 62 between the chord a'b' of Dq and t±, T2 are equal. By Claim IBTTl we 
have \9\ — 9 a \ < a and |#2 — &b\ < a > an d the claim follows. □ 



The connection between Euclidean stability and Q-stability. Let e pq be a ta-long Voronoi edge of the 
Euclidean diagram, for t > 9, and let u, v denote its endpoints. Let D u and D v denote the disks centered 
respectively at u, v, whose boundaries pass through p and q, and let D be a disk whose boundary passes 
through p and q, so that D C D u U D v and the angles between the tangents to D and to D u and D v at p (or 
at q) are at least ma each, where m > 4. (Recall that the angle between the tangents to D u and D v us at 
least ta > 9a.) 





c 
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Figure 26: The homothetic copy . 



Let c and p denote the center and radius of D, respectively. Note that c lies on e pq "somewhere in the 
middle", because of the angle condition assumed above. Let denote the homothetic copy of Q centered 
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at c and scaled by p, so Q c is fully contained in D and thus also in D U L)D V , implying that Q c is empty — it 
does not contain any point of P in its interior. (This scaling makes the unit circle Dq bounding Q coincide 
withD.) See Figure |26l 

Expand Q c about its center c until the first time it touches either p or q. Suppose, without loss of 
generality, that it touches p. Denote this placement of Q as Q c . Let £ p denote a supporting line of Q c at 
p. We claim that the angle between £ p and the tangent t p to D at p is at most a. Indeed, let £~ , £+ denote 

the tangents from p to Q c °\ By Claim IBTTl the angles that they form with the tangent t p to D at p are at 

most a each. As Qc is expanded to Q c , these tangents rotate towards each other, one clockwise and one 
counterclockwise so when they coincide (at Qq) the resulting supporting line £ p lies inside the double wedge 
between them. Since r p also lies inside this double wedge, and forms an angle of at most a with each of 
them, it follows that £ p must form an angle of at most a with r p , as claimed. 

Since the angle between the tangent t p to D at p and the tangent r p to D v at p is at least ma it follows 
that the angle between £ p and r p is at least (m — l)a. A similar argument shows that the angle between £ p 
and the tangent r p to D u at p is at least (m — l)a. 




Figure 27: The homothetic copy Q c . 



Now expand Q c by moving its center along the line passing through p and c, away from p, and scale it 
appropriately so that its boundary continues to pass through p, until it touches q too. Denote the center of 
the new placement as d , and the placement itself as Q c >. Let D c / be the corresponding homothetic copy of 
D centered at d and bounding Q c i. See Figure |27l 

We argue that Q c i is empty. By Claim [5T2l £ p is also a supporting line of Q c > at p. Refer to Figure [29] 
We denote by x p and y p the intersections of the supporting line £ p with dD c > and dD v , respectively. We 
denote by z the intersection of dD c > and dD v that lies on the same side of £ p as q. The angle Zpzx p is at 
most a since by Claim 15.11 the angle between £ p and the tangent to D c ' at x p is at most a. On the other 
hand the angle Zpzy p is at least (m — l)a since the angle between £ p and r p at p is at least (m — l)a. So 
it follows that the segment px p is fully contained in D v . Since the ray ~zp meets dD v (at p) before meeting 
dD c ', and the ray 'zx p ' meets dD c ' (at x p ) before meeting dD v , it follows that dD c ' and dD v intersect at a 
point on a ray between ~zp and ~zxp~. 

Let £ q denote a supporting line of Q c > at q. By Claim [531 the angles between pq and the lines l p , £ q differ 
by at most 2a. Since each of the angles between £ p and the two tangents r p and T p is at least (m — l)a, it 
follows that each of the angles between £ q and the two tangents and r q to D u and Z)„, respectively, at q, 
is at least (m — 3)a. 
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Figure 28: The homothetic copy Q c 




Figure 29: The segment px p is fully contained in D v . The circles dD dl dD v intersect at a point on a ray 
emanating from z between zp and z Xp . 



Refer now to Figure [30] We denote by z' the intersection of D c i and D v distinct from z, and we denote 
by x q , y q the intersections between £ q and D c >, D v , respectively. An argument analogous to the one given 
before shows that Zqz'x q < a while Zqz'y q > (m — 3)a. It follows that the segment qx q is fully contained 
in D v and we have an intersection between dD c i and dD v on a ray emanating from z' between the ray from 
z' to q and the ray from z' to x q . 

Our argument about the position of the intersections between D c > and D v implies that the entire section 
of dD c i between x p and x q is contained D v . Therefore the portion of Q c i to the right of the line through p 
and q (in the configuration depicted in the figures) is fully contained in D v . A symmetric argument shows 
that the portion of Q c > to the left of the line through p and q is fully contained in D u . Since D u U D v is 
empty we conclude that Q c t is empty 

The emptiness of Q c i implies that p and q are neighbors in the Q-Voronoi diagram, and that c' lies on 
their common Q-Voronoi edge ep q . 

We thus obtain the following theorem. 
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Theorem 5.4. Let P, a, and Q be as above. Then ( i) every 9a-stable edge of the Euclidean Delaunay 
triangulation is an a-stable edge ofDT®(P). (ii) Conversely, every 9a-stable edge of~DT®(P) is also an 
a-stable edge in the Euclidean norm. 

Note that parts (i) and (ii) are generalizations of Lemmas 13.41 and 13.51 respectively (with weaker con- 
stants). 

Proof. Part (i) follows directly from the preceding analysis. Indeed, let pq be a ia-stable Delaunay edge, for 
t > 9, whose Voronoi counterpart has endpoints u and v. Let Q c i be the homothetic placement of Q, with 
center c! , that touches p and q. We have shown that Q c > has empty interior if the ray p = pd lies between 
pu and pv and spans an angle of at least 4a with each of them. Assuming t > 9, such rays p form a cone of 
size (t — 8)a > a, which, in turn, gives the first part of the theorem. 

Part (ii) follows from part (i) by repeating, almost verbatim, the proof of Lemma 1331 □ 

There are many interesting open problems that arise here. One of the main problems is to extend the 
class of sets Q for which a near quadratic bound on the number of topological changes in DT Q (P), under 
algebraic motion of bounded degree of the points of P, can be established. 
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